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1  EXECUTIVE  SUMMARY 


The  overall  objective  of  this  three  year  research  project  is  to  develop  antenna  and  sensor 
systems  that  will  meet  the  Air  Force's  ever-increasing  demand  for  reduced  cost,  size,  weight,  and 
power  (CSWAP)  while  maintaining  or  increasing  system  functionality.  This  basic  research  effort 
consists  of  several  sub-tasks  that  use  analytical  and  numerical  methods  to  examine  the 
interaction  of  metamaterials  and  antennas,  topics  in  transformation  optics,  and  performance 
limitations  of  antennas  on  artificial  magnetic  grounds  and  unmanned  aircraft  systems  (UAS). 
Specifically,  this  in-house  research  and  development  program  consists  of  four  sub-tasks  which 
are:  1)  Metamaterial  Inspired  Superstrates  for  Smart  Antenna  Systems  (Dr.  Jeffery  Allen).  2) 
Bandwidth  Limitations  of  Antennas  on  Artificial  Magnetic  Grounds  (Dr.  Herscovici).  3)  Topics 
in  Transformation  Optics  and  Metamaterials  (Dr.  Bae-Ian  Wu).  4)  Fundamental  Performance 
Limitations  of  Antennas  Mounted  on  Electrically  Small  UAS  (Dr.  Kramer).  The  objective(s)  for 
each  sub-task  is  briefly  stated  below. 

1 .  Subtask  1 :  Metamaterial  Inspired  Superstrates  for  Smart  Antenna  Systems 

•  Create  simplified  methods  that  combine  material  and  physical  modeling 
approaches  to  understand  how  the  behavior  of  phased  array  antenna  systems  can 
be  modified  using  metamaterial  superstrates. 

2.  Subtask  2:  Bandwidth  Limitations  of  Antennas  on  Artificial  Magnetic  Grounds 

•  The  objective  is  to  determine  and  compare  bandwidth  limitations  of  antennas  on 
artificial  magnetic  grounds  (AMG)  and  traditional  perfectly  conducting  (PEC) 
ground  planes  as  a  function  of  the  separation  between  antenna  and  the  AMG/PEC 
plane.  The  goal  is  to  examine  the  trade-off  between  bandwidth  and  separation  for 
the  two  types  of  ground  planes  that  are  common  to  low-profile  antennas. 

3.  Subtask  3:  Topics  in  Transformation  Optics  and  Metamaterials 

•  In  the  first  year  of  research,  the  objectives  from  the  proposed  graphene  study  are 
to  yield  a  better  understanding  of  electromagnetics  in  graphene-based  structures 
and  provide  new  theoretical  tools  regarding  the  guidance  properties  of  graphene- 
based  transmission  line  for  metamaterial  devices  and  sensors  in  the  mm-wave  and 
optical  frequencies. 

•  For  the  transformation  optics/metamaterial  research,  the  objectives  would  be  to 
establish  a  working  code  based  on  analytic  technique  for  the  study  of  the  antenna 
problem.  The  configuration  for  such  hybrid  analytic  studies  will  be  of  canonical 
shapes,  including  planar,  cylindrical,  and  spherical  with  or  without  the  presence  of 
ground  plane.  Conformal  antenna  based  on  transformation  optics  would  be  the 
next  step.  At  the  same  time,  the  feasibility  of  using  spectral  transformation  and 
surface  current  based  media  inversion  is  expected  to  be  completed  by  the  first 
year  mark. 

4.  Subtask  4:  Fundamental  Performance  Limitations  of  Antennas  Mounted  on  Electrically 
Small  Unmanned  Aircraft  Systems  (UAS) 

•  The  goal  of  this  effort  is  to  establish  performance  limits  for  antennas  mounted  on 
structures  (e.g.  UAS,  Micro-UAS,  etc.)  with  sub-wavelength  dimensions. 
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The  remainder  of  this  report  discusses  the  objectives  in  more  detail  and  summarizes  the  results  of 
the  first  20  months  of  this  three  year  project. 
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2  Subtask  1 :  Metamaterial  Inspired  Superstrates  for  Smart  Antenna  Systems 

In  compact  phased  arrays,  only  a  small  number  of  antenna  elements  are  available.  This 
directly  determines  and  limits  the  number  of  degrees  of  freedom  that  are  available  for  adaptive 
beam/null-forming.  It  limits  the  number  of  nulls/beams  that  the  system  can  steer,  where  these 
null/beams  can  be  placed  and  their  quality  as  well  as  the  gain  of  the  system  [1],  Therefore,  for  a 
given  processing  algorithm,  system  performance  can  be  improved  either  by  increasing  the 
number  of  antennas  available,  increasing  how  efficiently  the  aperture  is  being  used,  or 
preferentially  moving  available  performance  (e.g.  ability  to  suppress  jammers)  to  critical  areas  of 
scan  space  thus  improving  performance  in  those  regions.  The  use  of  quasi-optical  elements  to 
enhance  the  performance  of  phased  array  anti-jam  and  beam  steering  systems  has  been  and 
continues  to  be  of  great  practical  interest.  Recent  developments  in  engineered  materials  suggest 
even  greater  design  possibilities  for  overcoming  system  limitations  than  were  previously 
available  with  standard  dielectric  materials.  Metamaterial-based  devices  provide  additional 
degrees  of  design  freedom  (e.g.  inhomogeneity,  anisotropy,  tenability,  etc.)  in  defining  the 
precisely  spatially  distributed  material  parameters  and  the  ability  to  adjust  and  optimize  them  so 
as  to  balance  the  various  tradeoffs  in  the  system. 


2.1  Objectives 

Our  objective  is  to  use  numerical  approximation  methods  to  reduce  the  complexity  of 
metamaterial  superstate  designs.  The  research  will  focus  on  developing  a  simplified 
methodology  to  solve  the  inverse  problem  which  uses  the  system  behavior  described  by  the  far- 
field  radiation  characteristics  to  calculate  the  material  parameters  that  result  in  the  required 
performance.  We  will  develop  this  method  in  the  context  of  metamaterial  superstrates  that  give 
the  desired  system  performance  in  antennas,  phased  array  antenna  and  scattering.  This 
investigation  also  aims  to  show  how  metamaterial  superstrates  can  be  designed  to  improve  far- 
field  characteristics,  scattering  problems  and  overall  system  performance  of  antenna  systems. 


2.2  Technical  Summary 
2.2.1  First  Year  Summary 

During  2012  this  subtask  consisted  of  setting  up,  configuring  and  coding  the  software 
problem  in  Matlab  and  computational  electromagnetic  solvers  (Comsol  [2]  and  HFSS  [3]).  In 
addition,  we  developed  analytical  and  full-wave  electromagnetic  simulation  models  for  the 
scattering  problem  to  show  the  efficacy  of  using  basis  functions  as  a  method  to  optimize  material 
properties  that  result  in  a  desired  effect  such  as  a  predefined  far-field  pattern  or  reduction  in 
scattering  cross-section.  Conducted  detailed  literature  search  to  help  devise  novel  models  to 
accurately  determine  and  modify  material  properties  of  given  electromagnetic  structure.  We  also 
setup  up  requisite  hardware  (high  performance  computational  computer  for  multi-threading)  in 
conjunction  with  configuration  of  computational  software  to  maximize  performance  and  speed  of 
aforementioned  simulations. 
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Our  initial  problem  has  concentrated  on  controlling  the  scattering  of  cylinders.  Using 
basis  functions  to  represent  the  spatial  distribution  of  material  properties,  we  started  exploring 
the  inverse  problem  of  minimizing  scattering  or  shaping  the  far-field  scattering  profile  to  a  pre¬ 
defined  profile.  This  task  began  in  fiscal  year  (FY)  2011,  and  is  scheduled  to  be  completed  in 
FY13.  Next  year  we  plan  to  continue  working  on  the  scattering  problem  and  also  extend  the 
model  to  the  radiation  problem.  Specifically,  we  will  modify  the  current  setup  to  include 
calculation  of  metrics  (e.g.  total  scattering  cross  section)  that  measure  the  performance  of  the 
simulations.  We  will  then  extend  the  simulations  to  three  dimensions. 


2.2.2  Second  Year  Summary 

During  2013  this  subtask  consisted  of  two  major  objectives:  demonstration  of  a  passive 
optimized  component  (lens)  and  inclusion  of  active  source  into  the  passive  transformed  space  for 
a  cohesive  combined  transformation.  We  use  the  source  transformations  approach  to  manipulate 
antenna  arrays  such  that  the  physical  location/shape  is  dictated  by  the  system  (e.g. 
wing/fuselage)  and  the  electromagnetic  performance  is  defined  by  the  application.  Numerical 
simulations  how  an  array  with  sources  arranged  in  a  circular  manner  can  be  surrounded  by  a 
transformation  optics  medium  so  that  it  can  be  made  to  receive  and  radiate  as  if  it  were  a  linear 
array  of  uniformly  spaced  elements  and  vice  versa.  We  have  finished  setting  up,  configuring  and 
coding  the  software  problem  in  Matlab  and  computational  electromagnetic  solvers  (Comsol)  for 
full-wave  electromagnetic  simulation  models  where  the  shape  of  individual  sources  was  changed 
using  Source  Transformation  Optics  (STO).  This  first  set  of  simulations  provided  solutions  to  the 
ideal  case  and  valuable  information  into  the  upper  bounds  defined  by  the  physics  of  the  problem. 
This  information  will  be  used  in  optimization  routines  to  develop  appropriate  approximate 
Transformation  Optics  (TO)  media  with  practical  parameters  that  can  be  fabricated  and  measured 
using  standard  testing  procedures.  Several  research  methods  exist  to  design  novel  engineered 
materials  such  as  metamaterials  to  control  the  radiation  patterns  of  antennas/arrays  and  also  the 
scattering  of  objects  illuminated  by  radiation  to  give  some  desired  output.  These  include  the 
following: 

-  Develop  theoretical  models  to  advance  understanding  of  how  device  and  material  parameters  as 
well  as  geometry  of  structures  affects  device  performance 

-  Design  simulation  studies  to  analyze  the  effect  of  constituent  material  design  parameters  on 
field  patterns,  aperture  efficiency,  resonances,  scattering  cross-section  and  impedance  of  a  device 

-  Develop  appropriate  figures-of-merit  to  analyze  the  efficacy  of  the  design  approach  as  well  as 
the  computational  efficiency  of  the  modeling  method 

The  challenges  are  that  (i)  the  problem  has  no  closed  form  solution  and  therefore 
optimization  using  full-wave  electromagnetic  simulations  required;  (ii)approach  can  be 
computationally  intensive  and  time  consuming;  and  (ii)  application  of  simple  optimization  using 
in-built  routines  in  the  FEM  software  package  are  not  efficient  and  also  do  not  produce  the 
required  results.  The  solution  is  to  design  a  method  that  gives  desired  behavior  and  also  produces 
practical  device  parameters  in  a  computationally  efficient  model.  AF  mission  depends  on 
increasingly  more  sophisticated  devices  with  unique  properties  not  available  in  natural  materials 
that  can  be  accessed  using  application-specific  engineered  materials  and  metamaterials.  This 
work  has  direct  implications  to  the  application  of  metamaterials  to  reduce  the  vulnerability  of 
systems  to  threats  such  as  jamming  and  anti-satellite  weapons  used  in  Precision  Navigation  and 
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Timing  in  GPS-Denied  Environments.  An  example  of  spatial  distribution  of  material  properties 
to  control  desired  device  characteristics  such  as  far-fields  and  scattering  parameters  is  shown 
below. 


Antenna  array 


v 


Optimized  metamaterial 


Figure  1:  Example  Setup  of  Superstrate  is  a  Transmission  System 

Superstates  can  be  designed  for  radiation  pattern  control  and  to  steer  a  beam  moving  the 
position  of  the  source  using  deterministic  material  parameter  distribution  similar  to  a  Luneburg 
lens.  They  can  also  be  used  for  scattering  parameters  control.  A  cylinder  illuminated  by  a  plane 
wave  can  be  coated  with  an  engineered  material  to  reduce  its  total  scattering  cross  section. 
Transformation  optics  (TO)  defines  the  forward  problem.  This  methodology  for  designing  media 
to  control  waves  using  accurately  designed  permittivity  and  permeability.  Coordinate  transforms 
can  be  used  in  form-invariant  Maxwell  equations.  Demonstrated  examples  of  TO  include  Cloaks, 
beam  shifters,  and  concentrators.  One  salient  example  of  TO  is  the  electromagnetic  cloak.  An 
object  is  surrounded  by  a  metamaterial  that  guides  incident  waves  to  original  trajectory.  To  the 

external  observer,  wavefront  appears  undisturbed. 


Cylindrical  Cloak 


Diamond  Cylinder  Cloak 


Flat  Cylinder  Cloak 


Ref  ^Allen.  J.,  N.  Kundtz  D.  Roberts.  S.  A  Cummer,  and  D  R.  Smith.  Source  transformations  using  superellipse  equations.  Applied  Physics  Letters,  v.  94  194101.2009. 


Figure  2:  Demonstrated  Examples  of  Transformation  Optics  Based  Cloaks 
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All  the  approaches  discussed  solve  the  forward  problem.  The  inverse  problem  is  more  interesting 
from  a  practical  standpoint  but  also  more  difficult.  The  complex  space  with  many  parameters 
makes  the  problem  intractable.  This  problem  gets  worse  in  3D  and  can  be  very  computationally 
intensive  and  can  result  in  unrealistic  material  parameters.  Optimization  provides  a  route  to 
simple  solutions  to  these  problems.  Some  results  using  optimization  are  discussed  below.One 
way  to  solve  the  inverse  problem  is  to  divide  the  solution  space  (material)  into  differential 
volumes  and  use  an  optimization  routine  to  obtain  a  desired  result  such  as  a  reducing  the 
scattering  cross  section  of  an  object.  This  can  create  very  large  solution  spaces  that  are  time 
consuming,  computationally  intensive  and  very  difficult  for  the  most  advanced  optimization 
routines  which  can  get  stuck  in  a  local  minima.  One  solution  is  to  reduce  the  number  of 
parameters  yet  keep  precision  of  the  solution  and  output  behavior  required  using  smart 
optimization.  The  goals  of  our  approach  are  to  (i)  Solve  the  inverse  problem  to  get  desired 
behavior  in  two  different  scenarios  as  follows.  Specifically  attempt  to  change  the  scattering 
cross-section,  and  reduce  scattering  cross-section.  Secondly,  change  far-field  pattern  of  a  given 
source  to  mimic  that  of  another  predetermined  source.  We  achieved  this  by  reducing  the 
parameters  in  the  solution  space  and  constraining  solutions  to  practical  material  parameters.  We 
used  basis  functions  to  represent  the  solution  space,  i.e.  spatial  material  properties.  We  started 
with  Zemike  polynomials  (basis  functions)  as  shown  below: 

2D  Zemike  functions  : 


ZTM 

[Znx[i\d] J  "  '  '  [coslm# 

where 


K(r)=T- 


i-iriM  -*! 


n  =  0.1,2 _  In -ml  even 


R  :  radial  polynomial,  n  =  0,1,2,...  :index  represents  degree  of  function  and,  m  =  -n  to  +n,  (n-m) 
even,  order.  Functions  are  orthogonal  over  unit  circle, 

JJ  3TM)  Z~(r,  9)rMff  -  — i-. 

Ovar  thamit  circle  \  ) 


where  l/(2n+l) :  normalization  constant. 

Orthogonality  to  express  any  function  f(r,0)  on  circle  as  a  sum  of  Zemike  functions 


/M)=X  Xcw^M) 

*2=0  m=-n 

2D  Zemike  analysis  can  be  extended  to  apply  to  3D  material  property  distribution 
3D  method:  Use  angular  variable  with  the  Zemike  radial  function 

where  ^nl  :  radial  Zemike  function  with  spherical  harmonics 
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Y;"(S4)  =  N'l"P;"{cos(0))e" 


N”  =  j-  OT?  ,  p”:  Legendre  functions 

\  4/T 


(l  +  m)  ' 

with  /<«,  /«  <  /  and  (n-1 )  must  be  even  also  apply  to  the  3D  Zernike 


Calculate  fields 


Calculate  criteria 


Check  criteria 


Applying  Zernike 

Optimization  chooses  amplitudes  for  the 

functions 

harmonics  (amn) 

w 

(Matlab) 

Call  CEM  and  simulate/compute  scattered 
fields 

(Call  Comsol  API  from  Matlab) 


Calculate  the  objective  function 

(Matlab) 


Figure  3:  Optimization  Process  Flow  Chart 
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Figure  4:  First  10  Zernike  Functions 


What  follows  are  a  few  practical  examples  of  how  such  a  optimization  process  using  Zernike 
bases  function  can  be  used.  First  we  start  off  with  the  goal  of  reducing  the  scattering  cross- 
section  of  given  object  (cylinder). The  scattered  electric  field  from  plane  wave  incident  at  350° 
without  metamaterial  shell  as  shown  below. 
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Figure  5:  Scattering  from  a  Cylindrical  PEC 


The  material  parameters  in  r  <  R2  are  mapped  into  Rj<r<R2  using  the  following  transformation 
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Can  be 

mathematically 
expressed  as 


r'=R1+r(R2-R1)/R2  </>'=</> 

Mapping  function  to  preserve 
polynomial  shape  but  accommodate 
scatterer  in  the  center 


Figure  6:  Mapping  Function/Transformation 

The  resulting  polynomial  Zemike  plot  is  shown  below  where  a  solid  single  mode  is  mapped  to 
create  a  central  discontinuity  using  the  mapping  function. 


Mapped/transformed 

Polynomial  single  basis  function  for  mode 


Figure  7:  Mapping  Function  as  Applied  to  a  Single  Zernike  Mode 


To  minimize  scattering  cross-section  and  total  scattered  field  we  define  the  objective  function, 
calculate  fields  using  full  wave  EM  models  and  optimize  material  parameters  using  Zemike 
moments.  Here  the  objective  function  to  be  minimized  is  the  total  scattered  field, 

IjF  {(hVt  ,  2W  n  m 

Scattering  cross  -  section,  uzcal  =  ^ — ,  Total  scattered  field, 

\E™\8)\  J= o<=o 

where  E5cat :  Scattered  field,  Einc  :  Incident  field,  9 :  incident  angle,  r,<p:  cylindrica  1  coordinate  s,0  <^<2^,0  <0  <2n 

In  this  example,  we  start  with  a  cylinder  (radius=0.125k;  s=3)  scattering  a  plane  wave  incident  at 
350°.  The  scattered  electric  fields  and  scattering  cross-section  are  calculated  as  a  baseline  before 
an  optimized  metamaterial  shell  is  placed  around  the  object  to  reduce  scattering. 
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Electric  field  ^ 


Scattering  cross-section 


angle' 


Figure  8:  Baseline  Scattering  from  a  Dielectric  Cylinder 


What  follows  are  two  solutions  obtained  using  the  process  described  above  showing  of  how  the 
scattering  from  the  dielectric  cylinder  can  be  reduced  using  a  thin  optimized  metamaterial  shell. 
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Figure  9:  Constant  Dielectric  Metamaterial  Shell  of  Negative  Permittivity 
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Figure  10:  Inhomogeneous  Dielectric  Metamaterial  Shell 

Next  we  optimize  a  structure  with  a  larger  shell  that  relies  on  guiding  the  waves  around  the 
scatterer.  Here  4  Zemike  moments  were  used  to  optimize  and  represent  the  spatial  distribution  of 

hr 

M4  =  W 

To  compare  the  results  with  prior  work  we  starting  with  the  equations  of  an  approximate  cloak: 
Optimize  sz  and  pr  while  holding  p(p  =1.  In  this  example,  we  start  with  a  cylinder  (radius=A/3; 
PEC)  scattering  a  plane  wave  incident  at  350°.  The  scattered  electric  fields  and  scattering  cross- 
section  are  calculated  as  a  baseline  before  optimized  metamaterial  shell  is  placed  around  the 
object  to  reduce  scattering. 


f  R1  ' 

fr-R.' I 

ft  — 

-  *"? 

V  r  ) 

Figure  11:  Scattering  Cross  Section  from  PEC  of  Radius=X/3 
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We  use  the  basis  function  optimization  to  design  a  shell  with  inner  radius  A/3  and  outer  radius 
4A/3  that  reduces  the  cross-section  of  a  scattering  cylinder.  The  resulting  material  parameters  are 
shown  below 


Figure  12:  The  Material  Parameters  Obtained  by  Optimization  for  the  Metamaterial  Shell 


Without  shell  Optimized  shell  Approximate  Cloak 


Figure  13:  Results  and  Comparisons  of  Scattering  for  Approximate  Cloak  and  Optimized 

Cloak 


2.2.3  Third  Year  Summary 

During  2014  this  subtask  consisted  of  two  major  objectives:  demonstration  of  a  passive 
optimized  component  (lens)  and  inclusion  of  active  source  into  the  passive  transformed  space  for 
a  cohesive  combined  transformation.  We  use  the  source  transformations  approach  to  manipulate 
antenna  arrays  such  that  the  physical  location/shape  is  dictated  by  the  system  (e.g. 
wing/fuselage)  and  the  electromagnetic  performance  is  defined  by  the  application.  Numerical 
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simulations  how  an  array  with  sources  arranged  in  a  circular  manner  can  be  surrounded  by  a 
transformation  optics  medium  so  that  it  can  be  made  to  receive  and  radiate  as  if  it  were  a  linear 
array  of  uniformly  spaced  elements  and  vice  versa.  We  have  finished  setting  up,  configuring  and 
coding  the  software  problem  in  Matlab  and  computational  electromagnetic  solvers  (Comsol)  for 
full-wave  electromagnetic  simulation  models  where  the  shape  of  individual  sources  was  changed 
using  Source  Transformation  Optics  (STO).  This  first  set  of  simulations  provided  solutions  to  the 
ideal  case  and  valuable  information  into  the  upper  bounds  defined  by  the  physics  of  the  problem. 
This  information  will  be  used  in  optimization  routines  to  develop  appropriate  approximate 
Transformation  Optics  (TO)  media  with  practical  parameters  that  can  be  fabricated  and  measured 
using  standard  testing  procedures. 

The  media  transformation  equations  for  the  source  transformation  of  sources  ar-ranged  in 
a  circular  manner  can  be  surrounded  by  a  transformation-optics  medium  so  that  it  can  be  made  to 
receive  and  radiate  as  if  it  were  a  linear  array  of  uniformly  spaced  elements  and  vice  versa  were 
derived.  However,  the  expressions  relating  the  original  coordinates  to  the  desired  coordinates  do 
not  have  inverse  expressions.  Hence,  the  summation  notation  representation  of  the  problem  does 
not  have  a  direct  solution.  We  use  properties  of  the  Jacobian  matrix  and  matrix  algebra  to 
formulate  the  expression  for  the  transformed  tensors.  This  approach  does  not  require  an  explicit 
expression  for  the  inverse  relationships  of  the  coordinates.  In  the  following  expressions  are  used 
to  transforming  a  uniform  array  of  line  sources  to  a  non-uniformly  spaced  circular  array  of  line 
sources. 

x(r',<p',z')  =  a(r')  cos  0' 
y(r\  4>',z')  =  b(r')  sintp' 
z(r',  0',  z')  =  z' 

where 

a(r')  =  p  ~~  ~jz  (r'  -  R±),  R2  <r'  <  R1 
K2  tix 

R2  --j  —  8  d 

b(r')  =  — - - —  (r'  —  /?x)  +  —  +  S,  R2  <  r'  <  R1 

R2-R1  1  2  2  1 

As  there  are  no  inverse  relationships  of  the  forms 

r'  =  r'(x,y,z ) 

0'  =  tp'(x,  y,  z) 

z’  =  z\x,  y,  z) 

the  summation  notation  and  explicit  formulations  for  the  tensor  elements  cannot  be  used.  Instead, 
we  rely  on  the  properties  of  the  Jacobian  matrix  A  where  A  is  defined  below 

r  _  _ A6qIA ^ 

^r(Pz  det/1  det4 

where  e0I  =  e08lJ  is  the  isotropic  and  homogeneous  material  tensor.  The  Jacobian  matrix  A  is 
explicitly  expressed  in  terms  of  the  derivatives  of  as: 

rdx  dx  dxi-1 
dr'  d<p'  dz' 

_  d(r',0',z')  _  ~  d(x,y, z)  -1  _  dy  dy_  dy_  _ 
d(x,y,z )  d(r',(f)',z ')  dr'  dtp'  dz' 

dz  dz  dz 
.dr’  dtp1  dz' . 


13 

Approved  for  public  release;  distribution  is  unlimited. 


the  permittivity  can  be  determined  from 


,  =  Ae^f_  =  B^e^CB^V 

6r<t>z  det4  det5_1 

One  advantage  of  using  this  approach  is  that  the  resulting  expression  will  be  in  terms  of  the 
primed  variables.  A  disadvantage  is  that  the  perspective  of  the  differential  field  behavior  due  to 
differential  material  changes  is  lost.  Also,  the  matrix  algebra  does  become  quite  tedious. 

The  Jacobi  matrix  A  is  first  calculated  from  B_1  where 


rdx 

dx 

dx  i 

dr' 

d<pr 

dz' 

B  = 

dy 

dy 

dy 

dr' 

dip' 

dz' 

dz 

dz 

dz 

dr' 

dtp' 

dz'  - 

and 


t>22 


da(i_)  cos  — a(r')  sin  0'  0 


dr' 


9i,<"7^sin  0'  b(r')cos0' 


dr' 


bn 

^21 

0 


hn  — 


R , 


•COS  0' 


bi2  — 


R2  £?i 

R2 

(r'  —  /?x)  sin0' 


R2  Ri 


R2-j-8 


b>2i  — 


Ro  —  Ri 


■sirup' 


r,  d  c 

R2  ~  2  ~  8  .  .  „  d 


_  _  (r'  —  R1)+  —  +  S 

R2-Ri  2 


COS  0' 


bi2 

^22 

0 


^33  —  1 

The  matrix  inverse  can  be  computed  using  Cramer's  Rule.  The  result  is  not  shown  here,  but  the 
resulting  elements  of  the  relative  permittivity  tensor  expressed  in  Cartesian  coordinates  are 


4R2(- R2n  )  y2+d2x2+4d8x2+482x2 
1  x \R 2~R1' 

t XX  — 


2(d+25)(I^=y-3/2v*2+y2 


I^x2(d+2V 


€ xy  €yx 


2R1y[  4x2y2+2x4+2y4-2R1y2 


Jx2+y2^j 


(  r2  r 


x3(d+2S)2- Jx2+y2 


..3/2 


yjx2+y2 


R2-Ri 


x2(d+2S)\ 


eyy 


4(  2x2y2+x4+y4-R1y2Jx2+y2^J 


(  R2  )2x2' 

x4  (d + 2  5) 2  •  (x2 +y  2 ) 

\R2-Ri)  x 

3/2 


yjx2+y2 
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^77 


2'3/2V^2+y2 


Lastly,  the  sources  are  also  transformed  using  the  Jacobian  matrix.  The  steps  are  not  shown 
here,  but  the  following  graphs  show  the  results  using  COMSOL  Multiphysics. 


JHUUkfl  p»rma*Oilty,  sas -i9mp$rv*nt  (1] 


Surf*ct:  ntlittv*  p*rmiaWrty.  My  cwwpoflint  til 


-S  -fi  -4  -2  0  2  0  B  T -2275,2 


SurffiGi:  [Mmflabilfty.  yy -(ompflrvfnt  (11 

o 


Figure  14:  Relative  Permeability  (same  as  relative  permittivity).  Due  to  Symmetry  of 
Tensor,  Only  Four  Elements  are  Shown.  Upper  Left:  j uxx,  exx ;  Upper  Right:  [iyy,  eyy; 
Lower  Left:  fixy  =  fiyx,  exy  =  eyx;  Lower  Right:  (izz,  ezz 
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Uniform  Linear  Array  in  Free  Space 
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Figure  15  Radiation  Patterns.  (Top)  Uniform  Linear  Array  in  Free  Space.  (Bottom) 

Circular  Array  with  Thick  Coating 

2.3  Conclusions  and  Future  Work 

In  summary,  we  specifically  addressed  the  following: 
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•  Scattering  inverse  problem:  Conducted  full  wave  simulations  using  Comsol  and  Matlab 
program  with  optimization  routines  to  include  basis  function  coefficients 

-  Small  dielectric  scatterer 

-  Approximate  cloak 

•  Radiation  inverse  problem  exploration 

-  Started  with  3  source  known  example 

-  Helped  research  idea  on  optimizing  material  parameters  using  basis  functions  to 
achieve  desired  far-field  pattern 

Next,  we  suggest  the  following  research  direction:  (i)  Derive  analytical  solutions  based  on  form- 
invariant  Maxwell’s  equations  in  the  STO  framework  for  the  transformed  sources  as  well  as  TO 
media;  (ii)  Develop  in-house  Matlab  code  to  simulate  the  passive  TO  case  operation  that  can 
provide  insight  into  the  wave  propagation  in  the  designed  media;  (iii)  Perform  finite  element 
method  (FEM)  based  simulations  using  software  packages  such  as  Comsol®  to  predict  the 
behavior  of  untransformed  and  transformed  arrays;  and  (iv)  Compare  results  from  (iii)  with 
respect  to  EM  fields  and  far-field  patterns  to  demonstrate  the  efficacy  of  the  approach. 
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3  Subtask  2:  Antennas  Operating  in  Close  Proximity  to  a  Textured  Surface 

In  this  work  a  new  method  to  control  the  interaction  between  antennas  and  nearby  objects  is 
proposed.  It  involves  the  use  of  a  thin  textured  surface  similar  to  a  high-impedance  layer  that  is 
wrapped  around  the  objects  that  are  located  in  close  proximity  of  the  antenna  (such  as  ground 
planes,  various  platforms,  etc.).  The  textured  surface  acts  as  a  “diffuser”  spreading  the  reflected 
fields  from  the  close  objects  and  thereby  reduces  the  power  otherwise  directed  in  specular 
reflections.  This  method  allows  reducing  the  profile  of  antennas  working  in  close  proximity  to 
various  objects. 

3.1  Introduction 

The  requirement  for  low  profile  antennas  that  are  conformal  to  large  and  finite  bodies  is 
not  new.  Typical  examples  are  antennas  mounted  on  an  aircraft  pod  (Figure  16a).  The  desire  is  to 
minimize  the  profile  of  these  pods,  which  would  bring  the  antennas  very  close  to  the  aircraft 
body  (Figure  16b).  The  aircraft  body  can  be  of  any  size,  shape  or  materials,  but  when  the  antenna 
is  brought  too  close  to  it,  the  performance  of  the  antenna  will,  in  most  cases,  deteriorate. 
Antennas  operating  in  the  close  proximity  of  finite  bodies  present  significant  deficiencies,  due  to 
the  strong  interaction  between  the  two.  Any  energy  emanating  from  the  antenna  undergoes  a 
reflection  from  the  object  present  in  close  proximity,  resulting  in  distortion  of  the  radiation 
pattern.  Edge  effects  from  the  finite  object  also  contribute  to  beam  degradation.  Squint  in  the 
main  beam,  increased  cross-polarization,  and  high-level  sidelobes  are  only  some  of  the  negative 
effects  of  this  interaction,  which  impose  constraints  on  the  design  of  a  low  profile  antenna 
system  that  almost  blends  with  the  finite  body. 


Figure  16  (a)  Traditional  High-Gain  Antenna  in  a  Pod;  (b)  Low  Profile  Antenna  in  a  Pod 
with  a.  Platform,  b.  Pedestal,  c.  High  Gain  Antenna,  d.  Radome,  e.  the  Direction  of  the 

Main  Beam 

This  work  describes  an  approach  to  address  the  distortion  of  the  pattern  due  to  the 
interactions  with  the  finite  body.  In  this  approach,  a  textured  periodic  surface  is  used  to  reduce 
the  interaction  such  that  the  negative  effects  are  largely  eliminated  over  a  certain  frequency  band. 
The  only  positive  effect  is  the  increased  directivity  resulting  from  the  array  of  the  antenna  with 
its  image. 


3.2  Problem  Definition 

The  proposed  test  case  is  shown  in  Figure  17.  An  end- fire  (low-profile)  dielectric  polyrod 
antenna  operates  in  free  space.  The  details  of  the  antenna  geometry  are  given  in  Table  1.  For 
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reference,  Figure  18  shows  its  radiation  pattern  in  the  elevation  (xy)  plane.  The  free  space 
radiation  pattern,  as  expected,  is  symmetric,  with  a  well-behaved  main  beam  and  a  gain  of  13.4 
dB  at  32GHz. 


Table  1  Antenna  Design  Geometry 


PEC  ground  size 

217x217mm  (23x23  X  @  32GHz) 

Antenna  Length 

29.25mm  (~3  X  @  32GHz) 

Distance  from  the  PEC  to  the  axis  of  the  antenna 

8.635mm  (0.92  X  @  32GHz) 

Figure  17  A  Waveguide  Fed  Multimode  Dielectric  Polyrod  Antenna 
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Figure  18  The  Radiation  Pattern  in  the  XY-Plane 


Next,  the  antenna  is  positioned  in  close  proximity  to  a  flat  PEC  object.  While  image 
theory  provides  some  insight  into  the  problem,  due  to  the  complexity  of  the  equivalent  sources 
(both  electric  and  magnetic)  required  to  adequately  represent  the  antenna  and  the  finiteness  of  the 
PEC  object,  the  image  theory  may  be  only  partially  be  used  to  understand  the  phenomenology. 
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Since  the  interaction  occurs  in  the  reactive  region  of  the  antenna,  none  of  the  other  analytical 
methods  are  applicable.  Therefore,  the  use  of  numerical  full-wave  methods  is  required  (in  this 
case  Microwave  Studio  by  CST  and  HFSS  by  Ansoft  were  used). 

The  results  of  the  Microwave  Studio  analysis  of  the  polyrod  antenna  are  shown  in  Figure 
19.  Figure  19  shows  a.  the  elevation  pattern  of  the  polyrod  located  at  0.927  above  an  infinite 
PEC  ground  plane  and  b.  the  elevation  pattern  at  the  same  height  above  a  rectangular  finite  PEC 
ground  plane.  As  expected,  when  placed  above  an  infinite  ground  plane,  the  maximum  of  the 
main  lobe  is  in  the  endfire  direction  (x-axis).  The  differences  between  the  finite  PEC  ground  and 
the  infinite  PEC  ground  are  evident:  in  addition  to  the  nulls  shown  for  the  infinite  PEC  case,  the 
beam  is  squinted  by  about  17  degrees  (~50%  of  the  42  degree  beamwidth)  and  many  additional 
ripples  are  noticed  in  the  pattern.  Although  the  finite  ground  is  237  by  237,  (which  can  be 
considered  an  electrically  large  ground),  the  edge  effects  are  still  significant. 


Comparison  between  Finite  PEC  GND  and  Infinite  PEC  GND  ( 3 2 G H z } 


Figure  19  A  Comparison  Between  the  Elevation  Radiating  Patterns  When  the  Polyrod  is 
Operating  in  the  Proximity  of  an  Infinite  PEC  Ground  and  a  Finite  PEC  Ground 


3.3  The  New  Concept  for  Pattern  Reconstruction 

A  number  of  techniques  have  been  used  in  the  past  to  reduce  these  unwanted  effects: 

1.  Using  absorbing  materials  on  the  edges  of  the  finite  objects.  This  significantly  affects  the 
efficiency  of  the  radiator. 

2.  Using  serrations  on  the  edges  of  the  finite  objects.  This  significantly  increases  the  size  and 
complexity  of  the  antenna  system. 

3.  Shaping  the  edges  of  the  finite  objects.  This  method  is  applicable  only  for  electrically 
large  bodies. 

The  method  proposed  here  attempts  to  diffuse  the  specular  reflections  from  the  scattterer. 
The  effect  of  diffusion  is  achieved  by  mode  conversion  of  the  incident  wave.  The  finite  object  is 
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covered  by  a  textured  surface  that  has  the  capability  to  reflect  and  convert  an  incident  plane  wave 
into  a  combination  of  higher  order  modes.  Specifically,  assuming  that  a  textured  surface  is  a 
periodic  structure  that  operates  in  the  grating  lobes  regime,  the  characteristics  of  the  unit  cell  are 
such  that  for  any  incident  plane  wave,  less  power  is  delivered  to  the  dominant  mode,  and  more 
power  is  delivered  to  the  higher  order  modes.  This  corresponds  to  a  lower  level  specular 
reflection  pattern  and  a  higher  envelope  of  minor  reflection  sidelobes. 

For  example,  Figure  20  shows  a  textured  surface  unit  cell  called  higher  order  mode  diffuser. 
Note.  The  geometry  looks  similar  to  a  typical  EBG  unit  cell,  however  the  functionality  is 
completely  different.  The  only  thing  in  common  between  the  EBG  unit  cells  and  this  higher 
order  mode  diffuser  unit  cell  is  the  fact  that  they  are  both  periodic  structures.  The  dimensions, 
the  shape  and  the  frequency  band  of  operations  are  such  that  the  unit  cell,  for  any  incident  plane 
wave  it  supports  and  excite  higher  order  modes.  Figure  21  shows  the  electric  field  distributions 
for  some  of  the  higher  order  modes. 


Figure  20  The  Higher  Order  Mode  Diffuser  Unit  Cell 
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Figure  21  The  Field  Distributions  of  Various  Higher  Order  Modes 
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The  analysis  of  this  particular  unit  cell  shows  that  for  a  normal  incidence,  the  power  delivered  to 
mode  0  (specular  reflection)  is  somewhat  reduced  and  is  comparable  to  the  level  of  power 
delivered  to  mode  6  and  7  (Figure  22).  This  means  that  the  power  delivered  to  the  main 
reflection  beam  is  reduced.  The  effect  is  significant  at  32  GHz,  where  the  power  delivered  to 
mode  6  and  7  is  comparable  to  the  power  delivered  to  mode  0.  The  effect  is  negligible  at  40 
GHz,  there  even  thought  ‘some’  power  is  delivered  to  mode  6  and  7,  the  level  is  not  high  enough 
to  produce  the  desired  effect.  As  shown  in  Figure  23,  the  desired  effect  is  achieved  in  three 
frequency  bands;  in  all  cases  the  reduction  in  power  of  the  dominant  mode  is  higher  than  ldB. 
Considering  that  the  main  beam  destruction  noted  in  Figure  3  is  the  result  of  the  superposition  of 
a  direct  wave  and  a  reflected  wave,  the  attenuation  of  the  reflected  wave  by  more  than  a  1  dB  can 
avoid  the  destructive  superposition  between  the  direct  wave  and  a  reflected  wave.  This  way,  the 
detrimental  impact  of  the  finite  objects  is  considerably  reduced. 

Mode  Spectrum 


-♦-32GHz  -*-40GHz 


Mode  Number 


Figure  22  The  Power  Delivered  to  the  Various  Modes  at  32  GHz  and  40  GHz 
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This  principle  was  applied  to  the  waveguide  fed  multimode  dielectric  polyrod  antenna  located 
over  the  finite  ground  mentioned  above.  This  time,  the  ground  was  covered  by  an  array  of  unit 
cells  described  above  (Figure  24). 


lifiEESiFrrw- 


Figure  24  The  Dielectric  Polyrod  Located  at  0.92X  above  a  Finite  Textured  Surface  Backed 

by  a  PEC  Ground 

Figure  25  shows  the  elevation  radiation  patterns  of  this  combination  at  32GHz,  compared  to  the 
same  pattern  in  the  absence  of  the  textured  surface  and  the  pattern  of  the  element  operating  in 
free  space.  As  shown  in  Figure  25,  in  the  presence  of  the  textured  surface  the  polyrod  antenna 
radiation  pattern  is  again,  well  behaved.  No  ripples  are  noticed  in  the  main  beam  and  the  sidelobe 
levels  are  close  to  the  ones  for  the  free  space  case.  The  presence  of  the  ground  is  noted  by  the 
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gain  in  directivity  as  predicted  by  the  image  theory.  As  mentioned  before,  the  squint  is  due  to 
the  finiteness  of  the  textured  surface  and  is  not  addressed  at  the  moment  by  this  technique. 
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Figure  25  Elevation  Patterns  for  the  Dielectric  Polyrod  in  Free  Space,  Located  at  0.92A, 
above  a  PEC  Ground  and  Located  at  0.92X  above  a  Finite  Textured  Surface  Backed  by 

PEC  Ground 


3.4  Conclusion 

A  new  method  for  improving  the  radiation  patterns  of  antennas  operating  in  close  proximity  of 
finite  objects  was  presented.  The  results  are  promising  and  further  investigation  is  required, 
mainly  regarding  the  bandwidth  capabilities  of  this  concept  and  the  robustness  of  the  method. 
The  inherent  edge  effects  might  be  reduced  using  a  non-uniform  textured  surface  with  specific 
characteristics  (such  as  surface  impedance  and  diffraction  coefficients),  which  may  require 
different  values  for  different  sections  of  the  surface.  This  is  a  topic  for  further  investigations. 
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4  Subtask  3:  Topics  in  Transformation  Optics  and  Metamaterials 

Metamaterials  and  transformation  optics  are  two  research  topics  that  are  of  significant 
importance  in  electromagnetics.  Transformation  optics  allows  us  to  find  the  required  constitutive 
parameters  associated  with  spatial  coordinate  transformation.  Metamaterials,  which  can  be 
described  as  synthesized  materials  with  constitutive  parameters  that  can  be  defined 
macroscopically  to  a  certain  practical  extent,  have  given  rise  to  a  branch  of  electromagnetics  in 
its  own  rights  and  have  provided  a  conceptual  means  for  realizing  the  complex  media  resulted 
from  TO. 


4.1  Objective 

In  recent  years,  transformation  electromagnetics  has  found  potential  applications  in  propagation, 
waveguiding,  scattering,  and  radiation.  For  antenna  applications,  using  transformation 
techniques,  one  can  transform  bulky  antennas  to  low  profile  ones.  In  general,  the  resulting 
medium  will  be  both  inhomogeneous  and  anisotropic.  In  this  section,  we  proposed  a  spherical 
core-shell  structure  which  can  achieve  arbitrarily  large  directivity.  We  investigated  the  problem 
by  finding  the  transformed  constitutive  tensors  and  solving  the  equivalent  problem  in  the  core¬ 
shell  configuration.  Using  the  Ricatti-Bessel  functions,  we  can  represent  the  field  components 
with  Debye  potentials  and  subsequently  solve  for  the  fields  in  all  regions.  We  applied  the 
formulation  to  several  cases  of  dipole  arrays  within  the  shell,  corresponding  to  both  free-space 
and  half-space  problems  in  the  virtual  space.  Overall,  the  calculation  demonstrated  that  the 
formation  of  virtual  aperture  is  indeed  theoretically  possible  and  the  effects  of  loss  on  the 
number  of  available  spherical  harmonics  and  directivity  are  investigated. 


4.2  Technical  Summary 
4.2.1  Introduction 

Transformation  electromagnetics  has  received  a  lot  of  interests  in  recent  years.  It  has  found 
potential  applications  in  several  major  areas  of  electromagnetics:  propagation,  waveguiding, 
scattering,  and  radiation.  Using  transformation  electromagnetics,  it  is  possible  to  manipulate  the 
electromagnetic  waves  to  achieve  novel  effects.  For  propagation  and  waveguiding,  transformed 
structures  allow  waves  to  propagate  and  bend  without  suffering  from  reflections.  Regarding 
scattering,  both  monostatic  and  bistatic  scattering  can  be  reduced.  In  this  section,  we  are 
primarily  interested  in  antenna  applications.  For  antenna  applications,  transformation  techniques, 
including  compression  and  conformal  mapping,  allow  the  transformation  of  bulky  antennas  to 
low  profile  ones.  For  example,  consider  the  parabolic  antenna  in  Figure  26  The  parabolic  antenna 
on  the  left  with  isotropic  medium  can  be  transformed  into  a  low  profile  antenna  on  the  right  [4], 
The  resulting  medium  will  be  both  inhomogeneous  and  anisotropic. 

For  the  horn  antenna  shown  in  Figure  27,  the  structure  on  the  left  with  isotropic  medium  can  be 
transformed  into  a  smaller  package  with  compression.  The  resulting  medium  in  this  case  is 
anisotropic  but  can  be  homogeneous.  This  has  also  found  potential  applications  for  decreasing 
the  spacing  between  antennas  and  ground  planes  to  reduce  the  typically  quarter-wavelength 
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distance.  The  resulting  medium,  compressed  in  the  normal  direction,  will  have  constitutive 
parameters:  s  =  xxaz\  +yyaz x  +z£(l/a)si  and  p  =  xxa\i\  +yya\x x  +zz(l/a)pi,  where  the  original 
medium  has  isotropic  permittivity  £1  and  permeability  p,.  The  thickness  is  reduced  by  a  factor  of 
a.  Notably,  the  idea  of  using  materials  with  high  permittivity  or  high  permeability  can  be  thought 
of  as  an  approximation  of  the  ideal  transformed  case.  Such  an  approximation  will  naturally  have 
good  performance  for  radiation  in  the  broadside  direction,  but  for  other  directions  equivalence  of 
the  complex  fields  in  the  entire  region  of  interest  with  respect  to  the  uncompressed  case  can 
only  be  obtained  via  the  inclusion  of  effective  anisotropy,  even  though  this  specific  transformation 
is  not  unique. 

Whether  we  are  transforming  a  parabolic  reflector  or  achieving  low  profile  via  simple 
compression,  the  effective  aperture  does  not  increase  without  space  projection.  This  can  be 
understood  from  the  previous  idea  of  “perfect  lens”  [5],  which  has  been  reinterpreted  using 
coordinate  transformations  [6].  Effectively,  space  projection  can  create  or  cancel  out  space  and 
maintain  the  perfect  matching  between  the  media. 


Figure  26  Transformation  of  a  Parabolic  Antenna  into  a  Low  Profile  Antenna.  Left: 
Original  Space.  Right:  Transformed  Space 


Figure  27  Transformation  of  a  Horn  Antenna.  Left:  Original  Space.  Right:  Transformed 

Space 


In  this  subtask,  we  present  a  layered  spherical  structure  that  is  able  to  form  a  virtual  aperture  in 
empty  space  by  applying  the  appropriate  space  projection.  Subsequently,  the  radiated  fields  will 
be  emitted  as  if  they  are  from  the  virtual  aperture.  Reducing  the  dimensions  of  radiators  is  always 
of  great  interest  in  the  antenna  engineering  community.  It  is  well  known  that  there  is  no 
mathematical  upper  limit  to  the  directive  gain  of  an  antenna  of  a  given  size  [7],  However,  there 
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still  lacks  a  general  design  capable  of  producing  arbitrarily  large  directivity.  With  our  proposed 
model,  any  arbitrarily  large  directivity  can  be  theoretically  achieved.  The  underlying  concept  has 
been  demonstrated  in  [8].  Here  we  formulate  the  problem  with  a  difference  in  terms  of  the 
potentials  used  to  represent  the  fields.  In  addition,  we  also  analyze  the  limit  of  the  achievable 
farfield  resolution. 

In  the  following  sections,  the  formulation  for  the  transformed  constitutive  parameters  and  the 
solution  to  the  radiation  problem  will  be  given  and  we  will  look  at  the  fields  and  farfield 
radiation  patterns  for  different  cases  and  explore  the  resolution  limit.  Whereas  many  problems 
in  transformation  electromagnetics  can  be  solved  numerically  via  electromagnetics  software 
packages,  it  is  hoped  that  the  analysis  here  will  provide  basis  for  future  comparisons,  allow  us  to 
analyze  the  effect  of  loss,  and  provide  some  insights  to  the  resolution  limit. 


4.2.2  Formulation 

In  this  section,  we  will  set  up  the  problem  in  such  a  way  that  instead  of  using  coordinate 
transformations  to  calculate  the  fields,  we  will  be  solving  a  transformed  version  of  the  original 
problem  in  the  physical  space.  The  transformation  will  result  in  a  shell  whose  medium  will  be 
radially  inhomogeneous  and  anisotropic.  We  will  be  using  the  time-harmonic  convention  e_l0J' 
throughout  this  section.  The  time  harmonic  Maxwell’s  equations  are: 

v  x~e=‘\(6b  , 

V  x  77=  -icoD+J  , 

V  D=p  , 

V  -77=0 , 

where  for  non  bianisotropic  media,  the  constitutive  relations  [9]  can  be  written  as 

D  =  £■  E  and  B  —  Ji  -  H  . 

The  permittivity  and  permeability  tensors,  in  general,  will  be  functions  of  space.  Figure  28 
shows  the  general  configuration  that  we  will  be  using  for  this  section.  Using  the  transformation  in 
[8],  the  empty  sphere  of  radius  c  is  compressed  into  a  smaller  sphere  of  radius  a  (region  1),  and 
the  empty  shell  between  r  =  b  and  r  =  c  is  mapped  into  another  shell  between  r  =  a  and  r  =  b  (region 
2).  Region  3  corresponds  to  the  space  where  r  >  b.  The  line  at  r  =  c  corresponds  to  the 
boundary  outside  which  is  the  region  of  equivalence.  To  solve  the  radiation  problem,  we  need 
to  find  the  transformed  constitutive  tensors  and  solve  the  transformed  problem. 


(1) 

(2) 

(3) 

(4) 

(5) 
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Figure  28  Transformation  of  Free  Space  into  a  Layered  Configuration.  Left:  Original 

Space.  Right:  Transformed  Space 


4.2.2.1  Transformation  of  the  Constitutive  Parameters 


Since  we  will  eventually  be  working  on  the  transformed  space  rather  than  the  original  virtual 
space,  in  this  subsection  we  will  use  the  unprimed  notation  to  signify  the  transformed  physical 
space  and  the  primed  notation  for  the  virtual  space.  For  other  parts  of  this  section,  the  primed 
notation  is  used  to  represent  the  source  coordinates.  Figure  29  shows  the  relationship  between  the 
radius  in  the  physical  space  r  and  the  radius  in  the  virtual  space  r1.  The  region  between  r  =  a  and  r 
=  b  has  a  negative  slope,  corresponding  to  space  projection.  The  resulting  constitutive  parameters 
are  thus  expected  to  have  negative  real  parts.  For  the  r  direction, 

r'  =  for  r  <  a  , 

r'  =  —  Ar  +  B  for  a  <  r  <  b.  and 

/  '  (o) 

r  —  r  for  r  >  b  . 


In  region  2, 


A  = 


c  —  b 


a 


and 


b(c  —  a 
b- 


(7) 


Under  this  transformation,  0  and  <p  are  unchanged,  thus  0=0'  and  cp  =  (p '. 
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r’ 


Figure  29  The  Relationship  Between  the  Radius  in  the  Physical  Space  r  and  the  Radius  in 

the  Virtual  Spacer’ 

Using  the  Jacobian  [10],  the  transformed  s  and  p  can  be  expressed  in  terms  of  the  original 
constitutive  tensors: 


_  =  ■ e  ■  Jn.,  and  _  =  Jn,  v  ■  Jn, 


\Jn< 


\Jrr' 


In  spherical  coordinates,  the  Jacobian  takes  the  form: 


dAr  1  dAr 

dr1  7  d0' 

rdAg  r  dAg 

dr1  7  dd' 

rsinddA^  /-sjn0  dA^ 

77  7~  ~7¥ 


1  dAr 
7 sin0'  d<j>' 

r  dAg 
r'sin©'  d(f)' 
r  dAq 

7  d<y 


(8) 


(9) 


The  resulting  Jacobian  is  a  constant  in  region  1  and  region  3,  and  it  is  diagonal  in  region  2: 


—  for  r  <a  , 

Jrri  =  diag[-^-,  jj,  jj)  for  a  <  r  <  b.  and 
J,,j  —  7  for  r  >  b  - 


(10) 


Using  these  results  in  (10),  we  have: 


£l  _  Mi  —  £ 

Gp  to  a 

6  =  |  =  diag[- 

fi  _  to  _  I 
Go  to 


— Ar+fl)- 
Ar- 


,  -A,  -A] 


for r  <a  . 

for  o  <  r  <  b.  and 
for  r>b  . 


(11) 


Note  that  in  the  above  expression  for  region  2,  A  and  B  are  positive  numbers  and  consequently 
the  constitutive  parameters  in  the  physical  space  for  the  shell  region  have  negative  real  parts. 
This  is  consistent  with  the  planar  case  in  which  space  projection  results  in  negative  permittivity 
and  permeability.  Because  we  will  be  analyzing  the  effects  of  loss,  we  introduce  a  factor  of  (1 
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-i8  )  to  both  s  and  p  in  region  2,  where  8  is  the  loss  tangent.  Note  that  in  region  1  both  8!  and  gi 
are  scalar  constants,  and  in  region  2  the  transverse  components  are  constant.  Only  the  radial 
components  in  region  2  are  functions  of  the  radius  r. 


4.2.2.2  Solution  to  the  Wave  Equation 

In  this  subsection  we  consider  a  radiator  inside  a  radially  inhomogeneous  and  anisotropic 
shell  as  a  result  of  the  coordinate  transformation  presented  in  the  previous  section.  This  problem 
can  be  solved  analytically  by  matching  the  boundary  conditions.  To  facilitate  this  effort,  we  will 
use  the  formulation  based  on  Mie  scattering  [9]  [10]  [11],  Consider  the  Hertzian  dipole  shown  in 
Figure  30,  oriented  in  the  a  direction,  and  positioned  at  (r1,  0 ',  <p ')  in  region  1 .  Decomposing  the 
waves  in  all  three  regions  into  TM  and  TE  components  with  respect  to  r,  we  can  express  the 
electric  and  magnetic  fields  in  terms  of  the  Debye  potentials  ne  and  nm.  For  TM  modes, 


Hjm  =  V  x  =  V  x  rne , 

£tm  =  — : —  (e  “ 1  •  V  X  V  X  rite)  • 

—lft) 


For  TE  modes, 

£TE  =  V  X  r<Pj£  =  v  X  77tm  , 

77tE  =  — (/d  _1  •  V  X  V  X  ?K, n)  . 

1  ft) 

We  can  also  use  the  potentials  to  represent  D  and  B  instead  of  E  and  H  [12]  [13]. 


(12) 

(13) 


(14) 

(15) 


Figure  30  Configuration  of  the  Problem.  The  Dipole  is  Oriented  Arbitrarily  in  the  d 

Direction  in  Region  1 
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Using  the  Maxwell  equations,  we  can  now  form  the  wave  equations  for  ne  and  nm.  Here,  we  will 
show  the  procedure  for  TM  modes,  noting  that  the  corresponding  results  for  TE  modes  can  be 
obtained  by  invoking  duality.  For  TM  modes,  the  wave  equation  can  be  written  as: 


vx— (r1 

-lft) 


•  V  x  V  x  rite)  =  icon  ■  V  x  rne  . 


(16) 


Note  that  only  the  transverse  components  of  the  permeability  tensor  will  be  used,  and  p  can 
effectively  be  treated  as  pt.  Also,  since  s  and  p  are  not  functions  of  r  and  are  constant,  we  can 
expand  both  sides  and  equate  the  arguments  of  the  curl  operators  with  the  addition  of  the 
gradient  of  a  scalar  potential  ®  to  obtain: 

i  dnte 


1 


0)2 


l  l _ !_ 

jj.,£r  r  r sin 


6 


1  1  0  t  drKf  \ 

P/£f  r  dA  30  > 

1  1 _ 1  d  /  dm,,  \ 

Ht£t  r  sind  dr  v  d<j)  ' 


d2rne 

d<j)2 


<■6 

+ 

1  d<Pe 

— 

r  d6 

1  d<P, 

rsinO  d(j) 

(17) 


For  0  and  cp  components, 

1  1  0  dme 

(O2  [i,£t  d6  dr 

We  can  thus  obtain  Cb  as: 


d<Pe 

do 


and 


1  1  d  d  rKe  0  <Pr 


<P,> 


1 


ft)2  f.lt£,  d<l>  dr 

1  dr  Tle 


0$ 


CO2  nt£r  dr 

Using  this  result,  we  have  the  following  for  the  r  component: 


£,  1  d2rne  1 


+ 


£t  r  dr2  r2  sin  6  d6 


d  .  .  d Tie. 
(sinO-T^r)  + 


1 


d2ne 


do 


r2sin2  6  d(j>- 


+  kf  —Tie  —  0  i 

£t 


(18) 


(19) 


(20) 


where  kt  is  the  wavenumber  corresponding  to  the  transverse  constitutive  parameters  with  respect 
to  r.  This  is  the  wave  equation  for  ne.  If  we  define: 


VP,=V>-1A(^) 


1 


r2  sin  6 


d  d%e  x 

5e(s,ne5s) 


+ 


d2Tle 


(21) 


r2sin2  0  d(f)2 


we  can  simplify  the  wave  equation  to: 

£,■  1  d2r7ie 


+  V2  ne  +  kj  —  lle  =  0  . 


£t  r  dr2 
In  particular,  when  zr=zt  and  [J,  =p,, 

(V2  +  k2)jt,.=  0, 

and  the  Helmholtz  equation  is  recovered.  This  is  the  case  for  region  1  and  3. 


(22) 


(23) 
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We  can  now  use  the  separation  of  variables  to  solve  for  ne.  We  do  so  by  letting  ne  =  /’(/')./o(0 14(  9 
).  Regarding /9,  we  have: 


d2/d 

d^2 


+  q2  f<p  =  0  and  /*  =  e±1,#  , 


(24) 


where  m  is  an  integer  from  0  to  n.  For  /e ,  we  have: 

1  d  /  •  d-fo  x  ,  ,  2 
(sinO— — )  +  (/r 


Q 


sin0d0v*  d  0 

In  terms  of  cos0 ,  with  p  =  n{n  + 1),  we  have: 

( 1  —  cos2  0 )  —  2  cos  0  -  ^'^6—  + 


)fe=  0. 


d(cos@): 


dcos# 


sin2  6 

n(n  +  1 )  — 


m 


1  —  cos2  6 


(25) 

fe  =  0, 

(26) 


for  which  the  associated  Legendre  polynomials  are  solutions  to.  Thus  /0  =Pm(cos0 ),  where  n  is  an 
integer  from  0  to  go.  Using  the  results  for  fp  and  /0,  regarding  /■  we  have:  n 

£r  d2rfr  [fr.2  2  /  |  i  \1  r  n 

+  —ktr  —  n{n  +  1 )  /r  =  0  . 

£t 


£r  1  dr2 


(27) 


At  this  point  we  can  put  in  the  explicit  r  dependence  : 


_  — 

~  r  —  R 


)2  ,  in  which  R  = 


b(c  —  a) 
0 c-b ) 


Rewriting,  we  have: 


d2rfr 
r~drr  + 


k2r2 _ - _ 

T  ( r-R )2 


«(«  +  1) 


fr  =  0  . 


(28) 


(29) 


If  for  fr  we  write  f  (r)  =  r~lF(k,r),  where  r  —  (r  -  R),  it  is  readily  shown  that  F  satisfies  the 
Riccati-Bessel  equation  [14]  [15]: 

(k,f)2F"  +  [(k,f)2  - n(n  +  l)]F  =  0 ,  (3#) 

where  F  =  b„  denotes  the  Riccati-Bessel  function,  defined  as: 

b n{ktr)  =  (ktr)bn(ktr)  ,  (31) 

in  which  b„  is  the  spherical  Bessel  function.  In  our  case,  we  only  use  the  spherical  Bessel 
function  of  the  first  kind  (j„)  and  the  spherical  Hankel  function  of  the  first  kind  (h(1)).  Because 
we  will  only  use  the  spherical  Hankel  function  of  the  first  kind,  we  suppress  the  superscript  and 
use  h„  from  here  onwards.  The  general  solution  for  ne/m  will  thus  be  a  double  summation  over  both 
m  and  n,  where  for  each  m,n: 


ne/m  =  -b„(k,r)P„  (cos 6 )e 


±i///0 


(32) 
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in  which  bn{k,r)  is  either  }„{k,r)  —  ( ktr)j„(ktr )  or  h „(k, r)  —  (k, r)h„(k, r) .  When  the  medium  is 

homogeneous  and  isotropic,  we  can  simply  set  k,=  k  and  R  =  0.  Throughout  this  section,  it  is 
useful  to  note  that: 


K 


j«(z)=i%Wz)’ 


K 


\Iy:  J„+Jr(Z)  +  iY«+^(Z) 


(33) 

(34) 


in  which  J  and  Y  are  the  cylindrical  Bessel  functions  of  the  first  and  second  kind  with  fractional 
order  n  +  1/2.  Also,  for  the  derivatives,  we  have: 

bln(z)  =  (z)bn-\(z)  -  nbn(z)  ,  (35) 


^p;"(cos0)  =  ^  [(/? cos 0 )PJ" (cos 6 )  —  (m  +  w)P„_ i (cos 0)]  . 


(36) 


These  expressions  are  useful  for  implementing  b„  and  evaluating  the  fields  numerically. 


4.2.2.3  Fields  and  Boundary  Conditions 


In  general,  the  fields  will  contain  both  TM  and  TE  components.  We  can  express  the 
electric  and  magnetic  fields  in  terms  of  ne  and  nm.  For  the  radial  components: 

„  1  ,d2rne  ,9  . 

Er  =  —  :  —  (  ..  ? — \-kt  rKe ) 


and 


Hr=7—( 


ic oer 

1  /d2r7l. 


m 


ico/J-r  dr 2 


+  kf  ritm )  • 


(37) 


(38) 


For  the  transverse  components,  using  the  notation  in  [1 1],  we  have: 

i  J_  d2  I  d 

C0£f  r  drdO  '  sind  d<p 
i  1  d~  „  d 

TV  n  W  '  V77 


~Ee~ 

E(p 

and 


'He 

H<!> 

1  d 

sind  dtp 
d 


l - 

bi 

1 _ 

Km 

i  1  d2 


(Wg,  r  drdO 

_ _ i _ 1  d2 

38  u>Hr  rsin0  drd<p ' 


7le 

K/n 

(39) 


(40) 


These  expressions  can  be  used  for  all  three  regions  to  find  the  fields  and  match  the  boundary 
conditions. 
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The  boundary  conditions  consist  of  matching  the  fields  at  the  source  and  the  interfaces.  For  the 
source,  which  in  our  case  exists  only  in  region  1,  we  use  the  expansion  of  a  Hertzian  dipole. 
Consider  a  dipole  moment.  For  the  TM  components,  we  have: 


— — a-V'x  V'xf' 

4k  k\ 


oo  n 

1 1 


n=Qm=-n 


(In  +  1  )(/7  —  |m|)! 
n(n  +  l)(n  +  |m|) ! 


(41) 


-j«(*r<)h»(*r> )Pjr‘ (cos a')?!;"1  (cos  9)em^~^  . 


For  the  TE  components,  we  have: 


n 


(2/7  +  1)(«  —  |  m  | ) ! 

„  7i(/7+  l)(n  +  \m\)\ 
(kr< )h„  (kr> )Pj7"'1  (cos  0')p»n|  (cos  0)eim^“^)  . 


H=0/n= 


(42) 


In  these  two  expressions  r’,  0’,  and  cp’  correspond  to  the  coordinates  of  the  source  and  the 
operators  V’  x  V’x  and  V’x  operate  only  on  the  source  coordinates.  Also,  r<=  min{ri,  r}  and  r>  = 
max{p,  r},  depending  on  the  radius  r  of  the  observation  point.  We  note  that  the  factor  P  is 
common  for  each  m  and  n.  Upon  suppressing,  we  can  simplify  the  notation  for  the  potentials 
corresponding  to  the  incident  spectrum  of  waves  in  region  1  as: 

(43) 

(44) 


•  /V 

K'\<  —  b\]n{k\r)  forr<r,and 


71 


i> 


A  . 

=ci\hn(k\r)  for  r>r  . 


Furthermore,  the  forms  for  TT™  and  ttte  are  identical,  therefore  we  can  simply  use  TTij2,3  for  each 
region,  keeping  in  mind  that  each  ny  2, 3  represents  two  cases.  Thus  for  region  1,  taking  into 
account  the  reflections  at  the  boundaries  [1 1],  we  write: 


n\<=a\R\2)n{k\r)  +  b\]n(k\r) 

for  r  <  r  , 

(45) 

K\>=a\  [hw (k\  r )  +  R i2j« (k\ r)] 

for  r>r. 

(46) 

In  region  2,  we  have: 

7t2  =  Cl2[hn{k2tr)  +  R23)n(k2tr)]  5 

(47) 

and  in  region  3, 

/■w  /X 

7t3  =  T\3Ci\K{hr) 

• 

(48) 

In  the  above  expressions 

involving  7ii>2j3, 

Tn 

(49) 

°2  \~R2\R2ia'  ’ 

p  _ rt  -  T2lR2lT\2 

l-RuRn 

7 

(50) 

and 
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(51) 


rfi  Tl2T23 

7 13  =  -  . 

1- *21*23 

Rn  and  f|3  are  the  total  reflection  and  transmission  coefficients,  which  account  for  the  multiple 
reflections.  The  reflection  and  transmission  coefficients  for  each  interface  can  be  found  by 
matching  the  boundary  conditions.  For  example: 

r7M 3/  /i  1  r  h/i  ( k2t  ci )  h„  (k  if  a )  —  ^/£\tjhthn(k2t  d)hn  ( k\tCi ) 

12  y/eu JU21  jn {k\ta)h'n (fo a)  -  y/£2tj^n{k2ta)}'n(kita) 

rTM  _  V^2tTUtjn{k2r^)in{kltCl)  —  \J ’£\tJ^)n(k\tCl))n(k2tCl) 

21  y/eitH2t]n(kita)f\'n(k2ta)  -  ^e^L^n{k2:a)]'n{k\ta)  ’ 


(52) 

(53) 


and 


tTM 
1 12 


i^2  t\J  M2r/£lr 

y/ Bit  Hit  'in  {kita)h'n  (kza)  -  y/e2tHit^n(,k2ta)yn{kita)  ’ 


tTM 

*2\ 


_ i£iry/ jUir/ffr _ 

y£]rM2/j n (k\ t a ) W„ (k2t a )  -  v^/jUih,, (fer fl )j„ (*i  / a ) 


(54) 


(55) 


We  will  also  need  the  coefficients  f?23  and  T23 ,  which  we  can  obtain  by  simply  changing 
subscripts  1  and  2  of  Rl2  and  Tn  to  2  and  3,  respectively.  To  obtain  the  coefficients  for  TE 
modes,  we  can  invoke  duality  by  swapping  s  and  p. 


We  now  derive  the  dipole  expansion  coefficients  for  a  dipole  corresponding  to  each  m,  n. 
Evaluating  the  coefficients  for  TM  and  TE  modes  yield: 


and 


where 


„™  =  =Amn(—0)  •  V'  x  V'  x  /s'b„(Jfei/)P!,M|(cose')e"iraf 


—A, 


k\  «cos0'a, 


r‘  sin  6 


i-bi(*,/)P!rl(cose')e-to*' 


+A„„  61(t,  /)p!r2 ,  (cos  e')e-‘”*' 


=  io)/i  1  Amn (  —  6)  ■  V'  x  ^„(ifcir')P!,m|(cos0/)e_i,nf 


1/77 


:  i(OHiAmn  -  -—prjbn(k\ r')P^ (cos  0')e 
r  sintT 


(56) 


(57) 
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(58) 


i  I£  (2n  +  l)(n—\m\)\ 

4 n  k\  n (n  -t- 1 ) (n  4-  \m\ ) ! 

We  have  now  completely  formulated  the  problem  for  the  potentials  and  fields.  In  the  next  section, 
we  will  apply  this  formulation  to  several  different  cases. 


4.2.3  Results  and  Discussions 


In  this  section,  we  will  show  the  computed  results  based  on  the  formulation  de-  scribed  in  the 
previous  section.  Subsequently,  these  computations  are  used  to  illustrate  various  aspects 
regarding  the  effects  of  enhanced  virtual  apertures.  In  particular,  we  are  interested  in  the  effects 
of  loss  on  resolution  in  the  farfield.  We  will  primarily  be  studying  the  0  -directed  electric  field 
in  the  xy  plane.  The  other  field  components  at  other  observation  points  can  be  obtained  in  a 
similar  fashion.  In  the  xy  plane,  0  =  n/2.  For  each  m  and  n,  the  electric  field  in  the  0  direction  is: 

Ee  =  +  ^(^f)pi-i(cos0)e''n<? 

U/Cf  / 

+«TEy  6i(fc#=)pJr,|(cos0)eim*  .  <*>> 

In  comparison,  the  corresponding  exact  expression  for  the  0  -directed  electric  field  in  the  xy  plane 
due  to  a  z-oriented  dipole  at  the  origin  in  free  space  is  [9] : 


Eq  = 


Qlkr 

4k  r 


1  ,  1  ,2 


l+Yr+{Tr) 


(60) 


Then,  based  on  (60),  we  can  calculate  the  reference  value  of  the  expected  electric  field 
corresponding  to  the  lossless  case. 


4.2.3. 1  4-Dipole  Array  in  Free  Space 

We  first  consider  the  configuration  as  shown  in  Figure  31.  In  the  corresponding  virtual  space, 
there  is  a  4-dipole  array  in  free  space.  The  dipoles  are  oriented  in  the  £  direction,  and  each  has  an 
in-phase  dipole  moment.  The  dipoles  are  separated  from  each  other  by  1/2  and  are  arranged  along 
the  x  axis.  In  the  transformed  space,  with  a  =  0.251 ,  b  =  0.51 ,  and  c=\l,  the  effective  physical 
aperture  is  reduced  from  c  to  b,  or  by  a  factor  of  2.  The  compression  ratio  ale  is  0.25 .  In  the 
ideal  case  with  no  loss,  we  expect  the  two  cases  to  produce  identical  fields  for  r>c. 


36 

Approved  for  public  release;  distribution  is  unlimited. 


Figure  31  A  4-Dipole  Array  in  Free  Space  with  X/2  Spacing.  In  the  Transformed  Space,  a 
=  0.25X,  b  =  0.5X,  and  c  =  IX,.  In  the  Lossless  Case,  the  Regions  for  r  >  c  are  Equivalent 


For  numerical  stability,  we  use  the  value  of  5  =  10  15  to  represent  the  lossless  case.  In  Figure  32 
we  show  the  field  value  \Re{E-  }  \  in  the  xy  plane.  The  range  is  between  ±6X  for  both  axes.  The 
main  beams  at  cp  =±n/2  as  well  as  the  sidelobes  are  clearly  observed.  In  Figure  33,  we  show  the 
field  value  \Re{Ez}  \  in  the  xy  plane  for  r>  b.  Visually,  they  look  similar  in  both  cases,  but,  as 
shown  in  Figure  34,  the  fields  are  in  fact  different  for  b  <  r  <  c.  Because  of  the  extremely  low- 
loss  and  the  well-matched  condition  created  by  the  transformation,  high  field  values  not  unlike 

that  of  surface  resonances  are  observed  at  points  near  r—b.  The  field  value  for  r>  c  converges 
to  the  reference  value  rapidly,  as  shown  in  Figure  35 . 


-6  -4  -2  0  2  4  6 

x(X) 

Figure  32  2D  Scan  of  |  Re{Ez}\  in  the  AT  Plane  Using  the  Exact  Expression.  The  Inner 
Circle  Corresponds  to  r  =  0.5X  and  the  Outer  Circle  Corresponds  to  r  =  IX 
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x(^) 

Figure  33  2D  Scan  of  |Re{Ez}|  in  the  XY  Plane  in  Region  3  (r  >  b)  Using  the 
Transformation  Method.  The  Circles  Correspond  to  r  =  0.25X,,  0.5^,  and  1  k 


Figure  34  Zoomed  2D  Scan  of  the  Magnitude  of  the  Difference  Between  the  Two  Methods. 
The  field  Values  for  r  >  c  are  Practically  Identical  While  the  Field  Values  for  b  <  r  <  c  are 
Different.  This  is  Expected  Because  the  Regions  of  Equivalence  are  for  r  >  c  only 
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Figure  35  Convergence  of  the  Electric  Field  Value  as  n  Increases 

In  Figure  36,  we  plot  the  magnitude  of  the  electric  field  component  Ee  at  r  =  100A,,  corresponding 
to  the  farfield  value.  The  pattern  is  effectively  identical  to  the  reference  one.  This  can  also  be 
seen  from  Figure  37,  where  we  plot  the  magnitude  of  the  coefficient  TTm  as  a  function  of  n.  This 
is  the  value  of  Tn  shown  in  section  2  for  TM  modes  given  n.  The  coefficients  Tte  have  the  same 
values  as  TTm  due  to  the  symmetry  between  82  and  g2-  In  the  lossless  case,  Ttm  =  0.25,  and  it 
remains  the  same  for  all  n,  suggesting  that  there  is  no  loss  of  angular  resolution,  and  the  original 
radiation  pattern  is  preserved.  Thus,  a  physical  aperture  with  diameter  X  has  formed  a  virtual 
aperture  with  diameter  2X  . 


4> 

Figure  36  Farfield  Radiation  Pattern  of  Ee  for  the  Physical  and  Virtual  4-Dipole  Array. 
The  Patterns  are  Essentially  Identical  Due  to  the  Low  Loss  Tangent  of  8  =  10-15 
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Figure  37  Magnitude  of  the  Coefficient  T  as  a  Function  of  n.  The  Magnitude  is  the  Same  as 
the  Compression  Ratio  a/c  =  0.25  for  n  from  1  to  20 


4.2.3. 2  4-Dipole  Array  with  Ground  Plane 

Consider  the  configuration  as  shown  in  Figure  38.  It  is  similar  to  the  one  in  Figure  31,  except 
that  here  the  dipole  elements  are  placed  X  /4  above  a  PEC  ground  plane.  In  the  transformed 
space,  a  =  0.25X,  b  =  0.5X,  and  c  =  IX.  The  compression  ratio  is  maintained  to  be  0.25.  Image 
theory  is  used  here  to  calculate  the  fields  where  we  put  four  dipoles  with  amplitudes  II  =  -1  at  y  = 
-0.0625X.  so  that  they  form  mirror  images  of  the  physical  dipoles.  The  region  of  interest  ranges 
only  from  (p  =0  to  n/2,  and  the  maximum  electric  field  strength  is  expected  to  be  double  that  of 
the  ffeespace  case.  In  Figure  39  we  plot  the  farfield  radiation  pattern  of  E{)  in  the  xy  plane.  The 
radiation  pattern  from  the  virtual  antenna  formed  by  a  smaller  core-shell  structure  is  equivalent 
to  that  from  a  physical  antenna  with  a  larger  dimension.  The  mainlobe  amplitude  and  the  sidelobe 
positions  agree  with  the  reference  case. 


Since  c  can  be  arbitrarily  large,  it  would  seem  that  we  can  theoretically  use  this  method  to  build 
a  structure  with  a  smaller  dimension  and  achieve  an  aperture  with  arbitrarily  large  directivity. 
However,  there  are  several  effects  we  need  to  consider,  including  coupling  and  loss.  By  virtue  of 
the  equivalence  due  to  the  transformation,  dipole  coupling  is  no  worse  than  the  situation  in  the 
virtual  space,  despite  the  spacing  between  the  elements  being  physically  much  closer.  However, 
loss  will  introduce  substantial  modifications  to  the  radiated  power  and  beamwidth.  This  effect  is 
shown  in  the  same  figure,  in  which  increasing  the  loss  causes  a  reduction  in  directivity  and  an 
expanded  beamwidth.  The  magnitude  of  the  beam  also  decreases  significantly  as  the  loss  tangent 
goes  up.  The  loss  also  introduces  both  mismatches  as  well  as  ohmic  loss  into  the  system. 
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In  Figure  40,  we  plot  the  magnitude  of  the  coefficient  T™  as  a  function  of  n  for  different  loss 
values.  It  is  observed  that  for  a  low  loss  tangent  of  5  =  1CT4,  the  transmission  coefficients  remain 
the  same  as  before  up  to  n  =  12.  As  the  loss  increases,  |T  TM|  decreases  more  rapidly.  For 
example,  when  5  =  10  2,  the  rolloff  occurs  at  n  =  6,  contributing  to  the  decrease  of  farfield 
resolution. 

In  Figure  41,  we  plot  the  farfield  radiation  pattern  of  Eq  in  the  xy  plane  for  c  =  2X  instead  of  IX. 
Hence,  the  compression  ratio  changes  from  0.25  to  0.125.  The  decay  of  performance  is  now  even 
more  rapid  than  in  the  previous  case.  For  the  same  5  as  before,  we  now  have  degraded 
beamwidth  and  magnitude.  This  is  also  evident  from  the  plot  of  the  transmission  coefficients  in 
Figure  42,  where  the  rolloff  occurs  for  smaller  n. 


y 


y 


Figure  38  A  4-Dipole  Array  with  X/2  Element  Spacing  Positioned  X/4  above  a  PEC 
Ground  Plane.  In  the  Transformed  Space,  a  =  0.25X,  b  =  0.5X,  and  c  =  IX.  Image  Theory  is 

Used  to  Calculate  the  Fields 
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Figure  39  Farfield  Radiation  Pattern  of  Ee  Corresponding  to  Various  8  for  the  Physical 
and  Virtual  4-Dipole  Array  Positioned  Over  a  Ground  Plane.  The  Compression  Ratio  is 

0.25  in  this  Case 


Figure  40  Magnitude  of  the  Coefficient  T™  as  a  Function  of  n.  The  Magnitude  Decreases 
as  n  Increases.  The  Decay  is  Faster  for  a  Higher  Loss  Value 
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Figure  41  Farfield  Radiation  Pattern  of  Ea  Corresponding  to  Various  5  for  the  Physical 
and  Virtual  4-Dipole  Array  Positioned  Over  a  Ground  Plane.  The  Compression  Ratio  is 

0.125  in  this  Case 
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Figure  42  Magnitude  of  the  Coefficient  T™  as  a  Function  of  n.  The  Magnitude  Decreases 
as  n  Increases.  The  Decay  is  Further  Increased  Due  to  the  Higher  Compression 


4.2.3.3  8-Dipole  Array  with  Ground  Plane 

Maintaining  the  same  compression  ratio  of  ale  =  0.125  and  element  spacing  in  the  virtual  space 
as  in  the  4-dipole  configuration,  we  now  change  the  number  of  dipoles  to  8,  keeping  the  same 
spacing  in  the  virtual  space.  In  Figure  43,  the  gray  dots  indicate  the  dipole  positions  in  the  virtual 
space  and  the  black  dots  correspond  to  the  transformed  ones.  When  we  had  4  dipoles  for  the  same 
compression  ratio,  we  were  not  taking  advantage  of  the  gain  in  virtual  aperture  size.  Here,  we 
are  now  taking  advantage  of  the  virtual  aperture  size  with  a  gain  factor  of  4,  from  a  physical 
aperture  of  IX  to  a  virtual  aperture  of  4X.  In  Figure  44,  it  is  seen  that  when  we  have  8  =  10“4,  the 
performance  is  discemibly  different  from  the  reference  case.  We  observe  a  more  dramatic  change 
in  the  radiation  patterns  for  each  corresponding  increase  in  8.  The  same  trend  for  | T  ™|  is 
illustrated  in  Figure  45. 
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Figure  43  An  8-Dipole  Array  Positioned  1/4  Above  a  PEC  Ground  Plane  with  7J2  Element 
Spacing.  In  the  Transformed  Space,  a  =  0.251,  b  =  0.5X,  and  c  =  2X.  Image  Theory  is  Used 

to  Calculate  the  Fields 


* 


Figure  44  Farfleld  Radiation  Pattern  of  £e  Corresponding  to  Various  8  for  the  Physical 
and  Virtual  8-Dipole  Array  Positioned  Over  a  Ground  Plane.  The  Compression  Ratio  is 

0.125  in  This  Case 
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Figure  45  Magnitude  of  the  Coefficient  T™  as  a  Function  of  n.  The  Magnitude  Decreases 
as  n  Increases.  The  Decay  Follows  the  Same  Trend  as  Before 

4.2.3.4  Virtual  Aperture  and  Farfield  Resolution 

Clearly,  there  is  a  relationship  between  the  transmission  coefficients,  the  beamwidth,  and  the  loss 
tangent  5.  In  this  subsection,  we  will  investigate  this  issue.  For  each  previous  case  where  the 
compression  ratio  is  0.125,  we  can  numerically  solve  for  the  null-to-null  beamwidth  of  Ee  by 
starting  from  the  normal  direction  y'  and  scanning  for  a  cp  that  corresponds  to  the  first  minimum 
field  value.  The  results  are  shown  in  Figure  46  and  Figure  47.  The  horizontal  axes  correspond  to 
the  loss  tangent  in  a  loglO  scale  and  the  vertical  axis  is  the  null-to-null  beamwidth  in  degrees.  In 
Figure  46,  for  the  4-dipole  case  subject  to  a  loss  tangent  smaller  than  10~3,  the  beamwidth 
converges  to  60°,  corresponding  to  the  lossless  case.  It  gradually  diverges  to  approximately  80° 
as  the  loss  tangent  reaches  10-1.  In  Figure  47,  for  the  8-dipole  case,  even  when  the  loss  tangent  is 
10-5,  the  beamwidth  does  not  converge  to  the  value  of  the  lossless  case.  It  also  gradually  diverges 
to  approximately  80°  as  the  loss  tangent  approaches  10-1.  However,  here  the  overall  angle 
deviation  is  larger.  For  the  same  range  in  8,  the  angle  now  changes  by  50°  rather  than  20°.  Even 
though  the  traces  of  the  plots  are  not  linear  with  respect  to  log10  8,  the  expansion  of  the  beam 
clearly  relates  more  closely  to  the  log  of  8  rather  than  just  8  itself. 

In  addition  to  observing  the  changes  in  beamwidth,  we  can  also  look  at  how  the  coefficients  T 
change  with  n.  T  are  the  TM/TE  transmission  coefficients  from  region  1  to  region  3  for  each  n. 
Their  complex  values  are  different  for  each  n  in  the  presence  of  loss.  In  Figure  48  and  Figure  49, 
we  show  the  3D  plots  of  T  as  a  function  of  n  and  8  ,  where  a  =  0.25X  ,  b  =  0.5X ,  and  c  =  IX  or  Tk. 
The  three  axes  correspond  to  n,  loglO  8,  and  T.  When  the  compression  ratio  is  0.25,  T  =  a/c  = 
0.25  for  the  lossless  case,  and  it  is  0.125  for  a/c  =  0.125. 

In  the  first  plot  (Figure  48),  with  a  compression  ratio  of  0.25,  the  flat  region  of  T  is  rather  large 
and  a  cliff-like  structure  is  observed.  The  edge  starts  approximately  from  n  =  8  on  the  n  axis  to 
loglO  8  =  -5  on  the  8  axis.  In  the  second  plot  (Figure  49),  with  a  compression  ratio  of  0.125,  the 
flat  region  of  T  is  smaller  and  a  cliff-like  structure  is  observed  again.  The  edge  starts 
approximately  from  n  =  8  on  the  n  axis  to  loglO  8  =  -10  on  the  8  axis.  These  two  plots  show 
quantitatively  how  the  transmission  coefficients  vary  with  n  and  8. 
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For  comparison,  consider  a  sphere  with  characteristic  dimension  kr,  where  k  is  the  wavenumber 
and  r  is  the  radius  of  the  sphere.  The  maximum  number  of  modes  supported  is  related  to  kr.  Take 
Figure  50  for  example,  where  we  are  plotting  the  same  parameters  as  in  Figure  48,  except  that  we 
use  an  isotropic  medium  for  region  2,  with  £2  =  £0  (1  +  id  ),  and  g2  =  po(l  +  id  ).  Effectively,  the 
sphere  has  a  radius  of  a  =  0.25k  despite  the  shell.  We  keep  the  medium  in  region  1  to  be  the 
same  as  the  previous  case  with  a/c  =  0.25;  thus,  ki  a  =  2n.  It  is  seen  that  T  drops  extremely  fast 
with  respect  to  n.  When  the  loss  tangent  is  high,  we  still  see  a  significant  drop  in  amplitude,  but 
we  no  longer  see  a  diagonal  cliff-like  shape.  This  is  one  reason  why  by  simply  putting  a 
radiating  source  inside  a  sphere  with  high  dielectric  medium  alone  will  not  support  high  angular 
resolution  beyond  a  certain  level  regulated  by  the  cutoffs  of  the  transmission  coefficients. 


Figure  46  Null-to-null  Beamwidth  for  a  4-Dipole  Array  with  a  =  0.25k,  b  =  0.5k,  and  c  = 
2k.  A  Smooth  Transition  is  Observed  as  5  Varies 


Figure  47  Null-to-null  Beamwidth  for  a  8-Dipole  Array  with  a  =  0.25k,  b  =  0.5k,  and  c  = 
2k.  The  Beamwidth  Broadens  Faster  than  the  4-Dipole  Case 
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Figure  49  3D  Plot  of  T  as  a  Function  of  n  and  5.  a  =  0.25X,  b  =  0.5X,,  and  c  =  Tk 


Figure  50  3D  Plot  of  T  as  a  Function  of  n  and  8.  a  =  0.25X,  b  =  0.5L  Instead  of  a 
Transformed  Shell,  a  Regular  Shell  from  r  =  a  to  b  is  used  with  zl  =  £0  (1  +  i5)  and  fi2  = 

ji0(1  +  i8) 
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Besides  studying  the  effects  of  the  transmission  coefficients  numerically,  we  can  also  look  at  it 


analytically.  Recall  that  for  each  n,  the  transmission  coefficient  Tu  can  be  written  as: 

rj,  _  TnTh 

13  ~~  i  -RuRx ' 


(61) 


In  the  above  equation,  Ti2,  T23,  R21  ,  and  R23  themselves  are  functions  of  the  constitutive  and 
geometrical  parameters.  Note  each  one  contains  the  Riccati-Bessel  functions  and  their 
derivatives.  We  can  look  at  the  arguments  of  the  Riccati-Bessel  functions  involved,  and  it  is 
straightforward  to  show  that  kia  =  koc.  Also,  k3b  =  kob.  Other  than  the  arguments  of  the  Riccati- 
Bessel  functions,  the  factors  S2tgi  and  83^2  are  the  same  as  Sig2t  and  S2t  p.3,  respectively. 
Therefore,  R21  =  R23  =  0  and  T12T23  =  a/c  for  the  lossless  case. 


Because  the  arguments  within  each  denominator  of  the  single  layer  coefficients  are  the  same  for 
the  lossless  case,  when  loss  is  present  the  denominators  of  R21  ,  R23 ,  Ti2  ,  T23  all  perturb  around 
their  lossless  Wronskian  value  of  “i”.  Whereas  the  nominators  of  T12  and  T23  remain  close  to  a/c 
under  perturbation,  the  nominators  of  R21  and  R23  are  extremely  sensitive  to  8.  Based  on  this,  it 
can  be  shown  that  within  the  scope  of  this  chapter,  the  contour  of  the  following  equation  yields  a 
rough  estimate  on  the  cutoffs  of  the  transmission  coefficients: 


|(k(«'i)k(«i)-k(«i)kW))  (h„(fi'2)hU«2)-h„(ff2)hU«2))|  =  1  .  (62) 

It  should  be  noted  that  for  planar  slabs  with  the  real  parts  of  the  constitutive  parameters  as  (-so, 
-go)  and  similar  loss  tangent  relationships  as  in  the  cases  shown  here,  the  spatial  resolution  in 
terms  of  the  maximum  tangential  wavenumber  before  transmission  cutoff  is  governed  by  a 
simple  expression:  (1/d)  ln(2/8  ),  where  8  is  the  loss  tangent  and  d  is  the  thickness  of  the  slab.  In 
both  the  planar  and  spherical  cases,  the  rapid  deterioration  of  the  spatial  and  angular  resolution  is 
related  to  the  log  of  8.  The  difference  is  that  in  the  virtual  aperture  case  considered  here,  the 
geometry  is  spherical  and  the  angular  harmonics  are  discrete  rather  than  continuous  like  the 
tangential  wavenumber  of  a  slab.  Consequently,  one  cannot  directly  solve  for  the  resolution. 
However,  one  can  plot  the  corresponding  contour  via  interpolation.  The  contour  plots  for  the  two 
virtual  aperture  cases  discussed  in  this  subsection  are  shown  in  Figure  51  and  Figure  52.  The 
contours  correspond  to  the  edges  of  the  cliff-like  shapes  in  the  previous  3D  plots.  It  is  seen  that 
the  contours  provide  reasonable  estimates  for  the  cutoffs  of  the  transmission  coefficients  and 
hence  the  farfield  resolution. 


n 


Figure  51  Interpolated  Contour  Plot  as  a  Function  of  n  and  8.  c  =  IX,.  The  Upper  Right 
Region  C  orresponds  to  the  Cutoffs  of  the  Transmission  Coefficients  and  Loss  of  Resolution 
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n 


Figure  52  Interpolated  Contour  Plot  as  a  Function  of  n  and  5.  c = 2k.  The  Upper  Right 
Region  Corresponds  to  the  Cutoffs  of  the  Transmission  Coefficients  and  Loss  of  Resolution 


4.2.4  Conclusion 

In  this  chapter,  we  proposed  a  spherical  core-shell  structure  which  can  achieve  arbitrarily  large 
directivity.  The  structure  was  predicted  and  obtained  from  the  proper  coordinate  transformations. 
A  large  virtual  aperture  can  be  projected  to  free  space  with  a  small  physical  dimension, 
subsequently  leading  to  a  large  directivity.  We  investigated  the  problem  by  finding  the 
transformed  constitutive  tensors  and  solving  the  equivalent  problem  in  the  core-shell 
configuration.  Using  the  Ricatti-Bessel  functions,  we  can  represent  the  field  components  with 
Debye  potentials  and  subsequently  solve  for  the  fields  in  all  regions.  We  applied  the  formulation 
to  several  cases  of  dipole  arrays  within  the  shell,  corresponding  to  both  free-space  and  half¬ 
space  problems  in  the  virtual  space.  Both  the  near  field  and  farfield  phenomena  were 
investigated.  The  relationship  between  the  beamwidths  and  transmission  coefficients  for  different 
loss  values  has  been  studied  and  an  estimation  for  the  cutoffs  of  the  transmission  coefficients  was 
provided.  Overall,  the  loss  is  linked  directly  to  the  farfield  resolution  in  terms  of  the  available 
angular  harmonics.  It  is  found  that  the  deterioration  of  beamwidth  changes  rapidly  with  the 
increase  of  loss  tangent  in  the  shell  medium.  Even  though  one  can  theoretically  construct  an 
antenna  aperture  with  arbitrarily  high  directivity  in  a  small  package,  based  on  the  specific 
transformation  used  here,  there  is  a  practical  limit  as  to  how  strongly  a  large  aperture  can  be 
com-  pressed  into  a  small  one.  It  can  be  reasoned  that  if  a  moderate  loss  tangent  value  can  be 
reached,  a  virtual  aperture  gain  factor  of  at  least  two  should  be  achievable. 
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5  Subtask  4:  An  Electric  and  Magnetic  Dipole  Antenna  for  Zero  Backscatter  and  Low 
Radiation  Quality  Factor 

The  original  goal  of  the  effort  was  to  establish  performance  limits  for  antennas  mounted  on 
structures  (e.g.  UAS,  Micro-UAS,  etc.)  with  sub-wavelength  dimensions.  However,  this  goal  was 
revised  because  other  researchers  achieved  a  significant  portion  of  the  proposed  research.  The 
new  goal  of  this  effort  is  to  use  analytical  and  numerical  methods  to  examine  the  possibility  of 
using  a  Huygens  source  (crossed  electric  and  magnetic  dipole)  to  achieve  a  low  Q  and/or  a  low 
backscatter  antenna. 


5.1  Objectives 

The  objective  is  to  use  analytical  and  numerical  methods  to  examine  the  possibility  of 
using  a  Huygens  source  to  realize  an  antenna  the  can  approach  the  Wheeler-Chu  Q  Limit  and  to 
realize  a  low  or  zero  backscatter  antenna.  While  the  main  emphasis  has  been  achieving  a  low 
backscatter  antenna  which  can  be  used  as  a  near  field  probe,  there  are  a  number  of  applications 
for  controlling  absorption  and  scattering  processes  from  antennas.  For  example,  a  Huygens 
source  can  also  be  used  to  realize  a  thin  layer  perfect  absorbing  material.  By  utilizing  antenna 
scattering  principles,  it  is  possible  to  design  absorbing  materials  at  high  frequencies  where 
magnetic  materials  are  not  possible  [16,17].  We  can  also  use  antenna  scattering  principles  in  the 
design  of  wireless  power  transfer  (WPT)  systems.  For  this  application,  antenna  scattering  can  be 
used  to  increase  power  transfer  which  is  important  consideration  in  wireless  charging  and  near¬ 
field  communication  (NFC)  devices  [18]. 


5.2  Technical  Summary 

Before  we  begin  to  discuss  the  design  and  realization  of  a  zero  backscatter  Huygens  source  we 
must  provide  an  overview  of  antenna  scattering  in  general.  In  the  following  sections  this  is 
accomplished  by  first  introducing  the  generalized  scattering  matrix  and  discussing  the  different 
ways  to  decompose  the  scattered  field  of  an  antenna.  We  then  proceed  to  discuss  the  basic 
approaches  for  controlling  scattering  and  the  important  role  that  absorbed  power  plays  in 
controlling  scattering.  After  concluding  the  overview  of  antenna  scattering  we  then  proceed  to 
discuss  scattering  from  an  ideal  Huygens  source.  This  is  followed  by  a  discussion  of  the 
realization  of  the  Huygens  source  and  several  implementation  issues. 

5.2.1  Generalized  Scattering  Matrix 

The  generalized  scattering  matrix,  S,  relates  the  outgoing  portion  of  the  total  field  to  the 
incoming  portion  of  the  total  field  by  the  following  equation 

rOutgoing\  _  $  flncoming\  (63) 

v  Modes  )  V  Modes  ) 

The  scattered  field  can  be  found  in  terms  of  a  S-matrix  by  noting  that  when  no  scatter  is 
present  the  scattering  matrix  reduces  to  the  identity  matrix  I.  In  this  case,  the  incoming  and 
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outgoing  modes  combine  to  give  the  standing  wave  representation  of  the  incident  field.  Thus  the 
scattered  field  of  any  object  is  given  by 

/ Scattered \  _  ^  (Incomings.  (64) 

V  Field  Modes  ) 

For  an  antenna,  the  generalized  scattering  matrix  is  slightly  more  complicated  than  a  typical 
target  because  all  antennas  have  a  port  which  connects  the  antenna  to  a  receiver/transmitter 
which  is  represented  by  some  load  impedance  ZL.  Therefore,  for  an  antenna,  the  generalized 
scattering  matrix  relates  incoming  and  outgoing  spherical  waves  (b;+  and  b,‘)  to  each  other  and  to 
the  signals  on  the  antenna  transmission  line  (a+  and  a')  which  is  illustrated  in  Figure  53.  In  Figure 
53,  Soo  is  the  reflection  coefficient  at  the  transmission  line-antenna  interface  looking  from  the 
termination  toward  the  antenna.  Similarly,  T  is  the  reflection  coefficient  at  the  transmission  line¬ 
load  interface  looking  from  the  antenna  toward  the  load. 


Figure  53  Illustration  of  the  Generalized  Scattering  Matrix  for  an  Antenna 

Based  on  the  illustration  of  Figure  53,  the  generalized  scattering  matrix  has  the  form  shown  in 
Figure  54.  From  Figure  54  it  is  observed  that  the  generalized  scattering  matrix  consists  of  the 
antenna  receive  modal  coefficients,  the  transmit  coefficients,  the  reflection  coefficient  and  the 
scattering  matrix.  Therefore,  the  generalized  scattering  matrix  completely  captures  the  radiation 
and  scattering  properties  of  an  antenna  which  makes  it  a  very  useful  tool  for  discussing  antenna 
scattering  for  reasons  that  will  become  apparent  in  the  following  section. 
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Figure  54  Generalized  Scattering  Matrix  for  an  Antenna 


5.2.2  Decomposition  of  Scattered  Field 

An  antenna’s  scattering  pattern  is  controlled  by  the  amplitude  and  phase  of  the  currents 
induced  just  like  any  other  structure.  However,  antennas  have  two  possible  modes  of  scatter. 
Firstly,  there  is  scattering  caused  by  the  shape,  size  and  material  (i.e.  structural  mode  or  residual 
mode)  and  this  mode  is  independent  of  how  the  antenna  receives  or  transmits  energy  (i.e.  doesn’t 
depend  on  transmit  or  receive  modal  coefficients).  Secondly,  there  is  scattering  based  on 
transmit/receive  properties  of  antenna  (i.e.  antenna  mode  or  load  dependent  mode)  and  this  mode 
depends  upon  mismatch  between  antenna  impedance  and  load  impedance  that  is  used  to 
terminate  antenna  port.  Therefore,  the  scattered  field  from  an  antenna  can  be  decomposed  into 
two  basic  terms  which  are  most  commonly  referred  to  as  the  structural  mode  and  the  antenna 
mode.  However,  we  will  refer  to  these  two  terms  as  the  residual  mode/term  and  the  load 
dependent  mode/term  because  the  commonly  used  terminology  often  leads  to  unnecessary 
confusion  in  the  author’s  experience. 

The  decomposition  of  the  scattered  field  can  be  written  in  terms  of  the  modal  components 
or  the  field  components.  For  example,  the  scattered  field  in  terms  of  modal  components  can  be 
written  as 


bs  =  bs°- 0.51(0X1 -rA)br  (65) 

In  equation  (65),  bs°  is  the  residual  mode  and  the  remaining  term  is  the  load  dependent  term 
which  depends  upon  the  short  circuit  current  1(0),  the  reflection  coefficient  rA  and  the  receive 
modal  coefficients  br.  Similarly,  the  scattered  field  can  be  expressed  in  terms  of  field 
components  as  follows 

ES(Z)  =  Es(Zl  =  0)  -  0.5/(0)(l  -  rA)Er  (66) 

In  equation  x,  the  first  term  on  the  right  hand  side  is  the  residual  mode  which  only  depends  on  Es 
and  corresponds  to  the  scattered  field  when  the  antenna  is  short  circuited.  The  second  term  on  the 
RHS  is  the  load  dependent  term  which  depends  upon  the  radiated  field  Er,  the  short  circuit 
current  and  the  reflection  coefficient. 
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An  important  aspect  concerning  the  decomposition  of  the  scattered  field  is  that  it  is  not 
unique  because  it  depends  upon  chosen  reference  impedance  that  is  used  to  terminate  the  antenna 
port.  There  are  three  standard  decompositions  that  are  commonly  found  in  the  literature 
[19,20,21],  They  are  the  short  circuit  (equation  (67)),  matched  load  (equation  (68))  and  complex 
conjugate  match  (equation  (69))  decompositions  and  are  given  by  the  following  equations. 

The  short  circuit  decomposition  (ZL=0)  is  given  by 

£S(Z)  =  Es{Zh  =  0)  -  0.5/(0)(l  -  rA)Er  (67) 

where  the  reflection  coefficient  is  defined  as 


71  zA  +  zL 

The  matched  load  decomposition  (ZL=ZA)  is  given  by 

ES(Z)  =  Es(Zl  =  ZA)  +  0.5/ (ZA)rAEr 

where  the  reflection  coefficient  is  the  same  as  the  short  circuit  decomposition. 
The  complex  conjugate  match  (ZL=ZA*)  decomposition  is  given  by 

ES(Z)  =  Es(Zl  =  Z/)  +  I(Z/)r/Er 

where  the  reflection  coefficient  is  define  as 


(68) 


(69) 


*  ZjA  —  Zj 

P  *  def  Zi _ 

A  ~  ZA+Z 

It  is  noted  that  throughout  the  remainder  of  this  report  that  the  complex  conjugate  matched 
decomposition  (equation  (69))  will  be  used  because  it  is  the  most  convenient  decomposition  to 
use  when  discussing  antenna  scattering  and  absorption. 

5.2.3  Controlling  Scattering 

5.2.3. 1  Two  Basic  Approaches  for  Controlling  Scattering 

One  of  the  main  advantages  of  decomposing  the  scattered  field  as  discussed  in  section  5.2.2  is 
that  it  provides  some  physical  insight  into  controlling  scattering.  All  of  the  decompositions  have 
two  basic  terms  which  are,  there  are  two  fundamental  approaches  for  controlling  scattering  (i.e. 
achieving  zero  backscatter)  [4,6].  They  are  the  matched  approach  (MA)  and  the 
compensation/cancellation  approach  (CCA). 

For  the  matched  approach,  scattering  is  minimized  by  minimizing  the  reflection 
coefficient  over  as  large  of  bandwidth  as  possible  which,  in  turn,  minimizes  the  load  dependent 
term.  The  residual  term  is  then  minimized  by  proper  design  of  the  antenna.  Ideally,  the  reflection 
coefficient  would  be  zero  over  the  desired  bandwidth  but  this  physically  impossible  to  achieve 
for  passive  lossless  matching  networks.  Furthermore,  this  is  not  achievable  with  active  matching 
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networks  because  of  stability  (i.e.  attaching  the  same  load  to  both  ports  of  a  non-Foster  matching 
network  results  in  an  oscillator). 

For  the  compensation/cancellation  approach,  scattering  is  minimized  by  using  the  load 
dependent  term  to  cancel  the  residual  term.  This  is  accomplished  with  proper  selection  of  the 
reflection  coefficient  which  can  easily  be  seen  by  rewriting  equation  (70)  as  follows 

ES(Z)  =  I(Z/)Er[Es(Z/)/I(Z/)Er  -  r/]  (70) 

We  can  then  define  the  first  term  inside  the  brackets  of  equation  (70)  as  A  which  results  in  the 
following  equation. 

ES(Z)  =  I(Z/)Er[A  -  T/]  (71) 

From  the  above  equation  it  is  apparent  that  the  scattered  field  can  be  minimized  by  minimizing 
the  expression  inside  the  brackets.  However,  the  extent  to  which  the  scattered  field  can  be 
controlled  with  this  approach  depends  upon  the  magnitude  of  A.  It  can  be  shown  [19]  that  if  the 
magnitude  of  A  is  greater  than  one  then  no  passive  load  can  reduce  total  scattered  field  to  zero, 
but  a  reactive  load  can  be  found  to  maximize  or  minimize  the  scattered  field.  If  the  magnitude  of 
A  is  equal  to  one,  then  there  is  a  passive  reactive  load  that  can  reduce  scattered  field  to  zero  [19]. 
When  the  magnitude  of  A  is  less  than  one,  there  is  a  dissipative  load  that  can  reduce  scattered 
field  to  zero  and  a  purely  reactive  load  that  will  maximize  scattering  [19]. 


5.23.2  Achieving  Zero  Backscatter  or  Zero  Total  Scattered  Power 

Both  approaches  outlined  in  section  5.2.3. 1  are  equally  valid  for  achieving  zero 
backscatter  or  zero  total  power  scattered.  For  the  matched  approach,  this  involves  matching  the 
antenna  such  that  rA*  =  0  and  designing  the  antenna  such  that  the  residual  mode  has  a  null  in  the 
backscatter  directions.  For  the  compensation/cancellation  approach,  zero  backscatter  can  be 
achieved  by  intentionally  mismatching  the  antenna  such  that  rA*  =  A.  However,  achieveing  zero 
back  scatter  or  zero  total  scattered  power  is  inherently  narrowband  regardless  of  the  approach 
used.  For  matched  approach,  a  reflectionless  matching  network  over  any  bandwidth  is  not 
physically  realizable  using  passive  and  lossless  devices  [22]  and  such  an  active  matching 
network  would  be  unstable.  For  the  compensation/cancellation  approach,  the  reflection 
coefficient  that  provides  minimum  or  zero  scattering  is  a  strong  function  of  frequency  and  aspect 
angle.  This  can  be  seen  from  equation  (70)  where  all  components  of  A  are  frequency  dependent 
and  the  field  components  are  angle  dependent.  This  makes  it  extremely  difficult,  if  not 
impossible,  to  achieve  zero  backscatter  over  any  finite  bandwidth.  However,  at  a  single 
frequency  it  is  possible  to  achieve  zero  backscatter  using  either  approach.  It  is  also  possible  to 
achieve  zero  total  scattered  power  but  this  can  only  be  accomplished  at  the  expense  of  absorbed 
power  which  we  will  discuss  in  the  following  section. 

5.2.3.3  Importance  of  Absorbed  Power 

Generally  speaking,  controlling  scattering  is  only  difficult  if  you  are  concerned  about  the 
absorbed  power,  in  particular  if  you  are  trying  to  simultaneously  maximize  the  absorbed  power 
will  minimizing  the  scattered  power.  The  reason  is  that  any  antenna  that  absorbs  power  must 
also  scatter  power  because  the  total  power  scattered  is  not  independent  of  the  absorbed  power. 
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They  are  related  to  one  another  through  the  forward  scattering  theorem  or  optical  theorem  [23], 
The  forward  scattering  theorem  (or  optical  theorem)  is  given  by  the  following  equation  which 
states  that  the  total  extinction  cross-section  atot,  which  is  equal  to  the  total  scattering  cross- 
section  as  plus  the  absorption  cross-section  ca  is  bounded  by  the  forward  scattering  amplitude 

f(0). 


An 

<*tot  =  os  +  oa  =  -ytmifiO)} 


(72) 


From  the  equation  above,  it  is  apparent  that  the  extinction  cross-section  cannot  be  zero  unless 
both  the  forward  scattering  amplitude  is  zero.  If  the  forward  scattering  amplitude  is  zero,  this 
means  that  the  scatter  or  antenna  does  not  interact  with  the  incident  field.  Therefore,  it  does  not 
scatter  the  incident  field  but  it  also  cannot  absorb  any  power  from  the  incident  field  because,  if  it 
did,  the  absorption  cross-section  would  be  greater  than  zero  which  leads  to  a  contradiction. 

There  are  two  basic  approaches  for  controlling  scattering  and  absorption:  1)  Minimize 
antenna  scattering  over  all  visible  space  (i.e.  minimize  total  scattered  power).  2)  Minimize 
antenna  scattering  over  a  limited  angular  region  of  space.  The  first  approach  involves  minimizing 
the  extinction  cross  section  which  is  subject  to  the  fundamental  bound  given  by  the  forward 
scattering  theorem  [23],  The  second  approach  typically  involves  reducing  the  bistatic  scattering 
from  the  antenna  in  the  back  half-space. 

For  the  first  approach,  minimizing  the  total  scattered  power  requires  the  antenna,  or  any 
object,  to  absorb  more  than  it  scatters  [19,24],  To  explain  this  fact  it  is  beneficial  to  make  use  of 
the  concept  of  absorption  efficiency  rj  which  is  define  as  follows 


V  = 


Oa  +  <*s 


(73) 


where  as  is  the  total  scattering  cross-section  and  aa  the  absorption  cross-section  which  are 
defined  as  follows. 


oa  =  An/Da 


(74) 


Os 


(75) 


In  the  above  equations,  Da  is  the  antenna  directivity  and  oi,  is  the  bistatic  scattering  cross  section. 
In  minimizing  the  total  scattered  power  we  are  effectively  maximizing  the  absorption  efficiency 
and,  ideally,  trying  to  make  rj  =  1.  In  order  for  the  absorption  efficiency  to  approach  one  we 
need  the  scattered  power  to  approach  zero.  However,  we  know  from  the  forward  scattering 
theorem  that  reducing  the  scattering  cross-section  also  reduces  the  absorption  cross-section 
because  we  are  decoupling  the  antenna  from  the  incident  field.  Therefore,  it  is  possible  for 
absorption  efficiency  to  approach  one  if  and  only  if  the  total  scattering  cross  section  approaches 
zero  faster  than  the  absorption  cross  section  [19,24],  This  results  in  a  trade-off  between  the  total 
scattered  power  and  total  absorbed  power.  That  is,  we  can  only  reduce  the  total  scattered  power 
(or  maximize  absorption  efficiency)  if  and  only  if  we  reduce  the  absorbed  power.  While  this 
approach  can  be  used  to  control  scattering  it  leads  to  antennas  that  do  not  absorb  as  much  power 
as  they  are  capable  of  which,  in  most  applications,  is  not  desirable.  However,  one  main  benefit  of 
this  approach  is  that  reducing  total  scattered  power  also  reduces  mutual  coupling  between 
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antennas.  Therefore,  this  approach  can  be  used  to  reduce  mutual  coupling  between  antennas  if 
desired. 
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Figure  55  Absorption  Efficiency  as  a  Function  of  Antenna  Directivity  and  Scatter 

Directivity  from  Reference  [24] 


For  the  second  approach,  we  are  interested  in  minimizing  antenna  scattering  over  only  a 
limited  angular  region  of  space.  This  makes  it  possible  to  reduce  scattering  over  a  certain  region 
without  reducing  the  absorbed  power.  Ideally,  the  power  scattered  in  the  back  half  space  would 
be  zero  while  the  absorbed  power  is  maximized.  For  antennas  that  are  capable  of  this  behavior 
the  total  scattered  power  is  greater  than  the  absorbed  power  and,  consequently,  the  antenna 
scatters  strongly  in  the  forward  direction.  In  practice,  scattering  can  be  reduced  over  a  limited 
solid  angle  with  minimal  to  no  reduction  in  absorbed  power.  This  method  relies  on  the  fact  that 
the  antenna  scattering  pattern  does  not  have  to  be  the  same  as  the  transmit/receive  pattern  which 
is  only  true  for  multi-mode  antennas.  This  concept  is  illustrated  in  Figure  56  which  shows  an 
antenna  whose  scattering  pattern  is  the  opposite  of  the  transmit/receive  pattern.  As  we  will 
discuss  in  the  following  sections,  a  Huygens  dipole  source  is  capable  of  this  behavior. 


Figure  56  Illustration  of  Controlling  Scattering  Using  Method  2 
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5.2.4  Generalized  Scattering  Matrix  for  an  Ideal  Huygens  Dipole  Source 


In  this  section  we  will  begin  to  discuss  scattering  from  an  ideal  Huygens  dipole  source  by 
first  deriving  the  generalized  scattering  matrix  for  an  ideal  Huygens  source.  This  will  be 
followed  by  a  discussion  of  the  scattering  properties  and  how  effectively  we  can  control  the 
scattering  from  the  Huygens  source  in  in  section  5.2.5. 

A  Huygens  dipole  source  consists  of  orthogonal  electric  and  magnetic  dipole  moments. 
This  type  of  antenna  is  one  of  the  simplest  multi-mode  antennas  and  it  is  possible  to  minimize 
scattering  base  on  method  1  (total  power  scattered/absorption  efficiency)  or  method  2  (limited 
solid  angle).  Also,  the  combination  of  electric  and  magnetic  dipole  moments  naturally  produces  a 
null  in  its  radiation  and  scattering  pattern  which  can  be  utilized  to  achieve  zero  backscatter. 

In  general,  only  for  simple  cases  (2  modes  or  less)  can  one  obtain  the  generalized 
scattering  matrix  analytically.  The  conditions  for  which  this  is  possible  are  there  can  only  be  two 
modes  per  port  in  the  generalized  scattering  matrix  and  the  symmetry  and  unitary  properties  must 
be  applicable.  Only  under  these  conditions  do  we  have  enough  degrees  of  freedom  to  determine 
the  generalized  scattering  matrix  analytically.  In  what  follows,  we  are  able  to  analytically  derive 
the  generalized  scattering  matrix  for  a  Huygens  source  because  we  assume  that  the  antenna  is 
reciprocal  and  lossless.  To  derive  the  generalized  scattering  matrix  for  a  Huygens  source,  let  us 
first  construct  an  ideal  Huygens  source  from  an  x-directed  electric  dipole  and  a  y-directed 
magnetic  dipole  as  shown  in  Figure  57. 
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Magnetic 
dipole 
moment 
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Electric  dipole 
moment 


Figure  57  Illustration  of  the  Geometry  for  a  Crossed  Electric  and  Magnetic  Dipole 

Moment 


We  can  construct  the  generalized  scattering  matrix  by  first  considering  the  transmit  case  where 
each  dipole  mode  is  excited.  We  consider  the  most  general  case  where  the  dipole  modes  have 
arbitrary  excitation  where  the  variable  k  is  the  ratio  of  the  excitation  or  modal  coefficients. 
Noting  that  there  is  a  90  degree  phase  difference  between  the  radiated  field  from  an  electric  and 
magnetic  dipole,  one  way  to  write  the  partial  scattering  matrix  is  as  follows 


/  0 

1 

V/c2  +  1 

Vk2  +  1 

-jk  | 

-jk 

i 

V  k2  +  1/ 

V  k2  +  1 

(76) 


where  k  is  the  ratio  of  the  modal  coefficients  and  is  real  (0  <  k  <  1)  and  the  modal  coefficients 
are  normalized  such  that  the  total  power  radiated  is  always  unity  which  is  the  origin  of  the 
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1/V/c2  +  1  term.  On  the  left  hand  side  (LHS)  of  equation  (76),  the  top  term  is  zero  because 
under  transmit  conditions  there  is  no  reflected  wave  on  the  transmission  line.  The  middle  term  in 
the  column  vector  represents  the  magnetic  dipole  coefficient  (mon)  and  the  third  term  the  electric 
dipole  coefficient  (n0n).  Using  the  symmetric  and  unitary  properties,  one  finds  that,  for  the 
matched  case,  the  S-matrix  must  be 
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13^13 
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^13 
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In  this  case  the  symmetry  property  ( Smn  =  Snm  for  m  =£  n)  is  a  consequence  of  reciprocity  and 
the  unitary  property  ([5] [5]*  =  [/])  is  a  consequence  of  energy  conservation  for  a  loss-less 
passive  system  (we  have  assumed  zero  ohmic  loss  in  this  case).  Substituting  the  S12  and  S13 
terms  from  equation  (76)  into  equation  (77)  we  obtain  the  generalized  scattering  matrix  for  an 
ideal  Huygens  source 
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In  equation  (78),  a  is  a  phase  angle  that  can  be  chosen  arbitrarily. 
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5.2.5  Scattering  from  an  Ideal  Huygens  Source 


In  this  section  we  will  obtain  the  scattered  field  from  an  ideal  Huygens  source  and 
discuss  the  minimization  of  the  total  scattered  power  and  the  necessary  conditions  for  zero 
backscatter.  We  begin  by  assuming  an  incident  field  of  the  form 


Et  =  xEie-iP* 


The  scattered  field  is  given  by 


Es  = 


C s-o 


(79) 


(80) 


By  substituting  equation  (77)  into  (80)  and  using  the  complex  conjugate  decomposition  of  the 
scattering  field  (equation  (69)),  the  residual  term  can  be  written  as 
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and  the  load  dependent  term  as 
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(82) 
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Now  that  we  have  the  scattered  field  from  and  ideal  Huygens  source  we  can  demonstrate  how 
scattering  can  be  controlled  using  the  methods  outlined  in  section  5. 2. 3. 3  for  the  cases  where  the 
antenna  load  is  the  complex  conjugate  of  the  antenna  impedance  and  the  case  where  it  is  not. 


5.2.5.1  Total  Scattered  Power  for  the  Complex  Conjugate  Matched  Case 

In  this  section  we  will  discuss  minimizing  the  total  scattered  power  when  the  load 
impedance  is  the  complex  conjugate  of  the  antenna  impedance.  We  want  to  minimize  total 
scattered  power  (max  absorption  efficiency)  for  a  load  impedance  equal  to  the  complex 
conjugate  of  the  antenna  impedance.  Therefore,  this  is  equivalent  to  simultaneously  minimizing 
the  total  scattering  power  while  maximizing  the  absorbed  power  or  maximizing  the  absorption 
efficiency.  For  a  complex  conjugate  match,  the  reflection  coefficient  is  zero  which  reduces 
equation  (69)  to  the  following 

E{Z)  =  E{Za *)  +  ri(Za*)Tr  ->  E(Z )  =  E(Za *)  (83) 

Therefore,  the  total  scattered  field  is  just  the  residual  term  from  equation  (69).  We  can  use  this  to 
then  compute  total  power  scattered  and  the  absorbed  power  as  follows.  The  scattered  power  is 
given  by 


PS  =  ^J&n~-frn+| 


(84) 


Using  equations  (78)  and  (81),  the  total  scattered  power 
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The  absorbed  power  is  given  by 
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(86) 


Pa  =  |a-|2 
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Using  equation  (78)  in  equation  (85),  the  absorbed  power  is  given  by 


^  3 CM)2  (k  +  l)2 
167771  k2  +  1 


(87) 


To  maximize  the  absorption  efficiency  we  note  that  scattered  power  must  be  less  than  the 
absorbed  power.  In  Figure  58  we  have  plotted  the  ratio  of  the  total  scattered  power  to  the 
absorbed  power  as  a  function  of  the  ratio  of  the  modal  coefficients  k  and  the  relative  phase  a. 
Because  Ps  must  be  less  than  Pa,  the  maximum  of  absorption  efficiency  must  be  located  in  the 
dark  blue  area  in  Figure  58.  It  turns  out  that  the  maximum  absorption  efficiency  occurs  for  k  ~ 
0.3  and  a  =  0.5ti  and  it  is  approximately  67%.  It  is  noted  that  the  point  where  the  absorption 
efficiency  is  maximized  corresponds  to  a  22.5%  reduction  in  absorbed  power.  That  is,  for  a 
Huygens  source,  maximum  absorbed  power  always  occurs  for  k=land  it  is  independent  of  a 
which  can  easily  be  seen  from  equation  (87).  It  is  possible  to  increase  the  absorption  efficiency  if 
the  antenna  is  not  complex  conjugate  matched  by  proper  choice  of  the  load  impedance.  Of 
course,  this  will  lead  to  a  further  reduction  in  absorbed  power  as  discussed  previously. 


k 


Figure  58  The  Ratio  of  the  Total  Scattered  Power  to  Total  Absorbed  Power  as  a  Function 
of  the  Ratio  of  the  Modal  Coefficients  k  and  the  Relative  Phase  a 


S.2.5.2  Zero  Backscatter  for  the  Complex  Conjugate  Matched  Case 

In  this  section  we  will  discuss  the  condition  for  zero  backscatter  when  the  antenna  is 
complex  conjugate  matched.  We  are  interested  in  the  zero  backscatter  condition  for  the  complex 
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conjugate  matched  case  because,  under  this  matching  condition,  we  have  maximum  absorbed 
power  or  power  transfer  to  the  load.  Again,  under  this  condition  the  reflection  coefficient  is  zero 
which  reduces  scattering  field  from  equation  (69)  to  the  following 

E(Z)  =  E(Za *)  +  n{Za*Wr  ->  E(Z)  =  E(Za*)  (88) 

The  residual  field  term  can  be  written  using  the  modal  functions  as  follows 

-  *  V3  EtX  (  ■  ■  c~  ,  c  \  (89) 

E(Za  )  =  — — j=  (0  m011  nell)l  J-JS22+S23 

4  vw  \-i-js32  +  s33J 

where  mon  is  the  magnetic  dipole  field  pattern  and  non  is  the  electric  dipole  field  pattern. 
Multiplying  the  column  vector  by  the  row  vector  in  equation  (89)  and  substituting  the  S- 
parameters  of  ideal  Huygens  source  we  have 
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Under  the  far  field  condition  (large  argument  approximation),  the  mode  functions  for  the 
backscattered  (0  =  n,  <p  =  0°)  field  are 
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Substituting  equations  (91)  and  (92)  into  equation  (90),  the  backscattered  field  is  given  by 
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From  equation  (93)  it  is  apparent  that  zero  backscatter  is  only  achievable  for  k=l  which 
corresponds  to  equal  excitation  of  the  electric  and  magnetic  dipole  modes.  For  this  particular 
case,  the  backscattered  field  is  not  affected  by  the  phasing  between  the  electric  and  magnetic 
dipoles.  This  result  is  similar  to  the  zero  backscatter  condition  for  small  isotropic  particles 
[25,26], 


5.2.5.3  Zero  Backscatter  for  an  Arbitrary  Load 

In  this  section  we  will  discuss  the  condition  for  zero  backscatter  when  the  antenna  is 
matched  to  an  arbitrary  load  impedance.  For  the  case  where  the  load  impedance  is  not  complex 
conjugate  matched  to  the  antenna  impedance  there  is  no  value  of  k  and/or  a  alone  that  can 
produce  zero  backscatter.  For  this  case,  zero  backscatter  can  only  be  achieved  by  proper 
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mismatch  which  leads  to  a  reduction  in  absorbed  power.  The  reflection  coefficient  for  which 
zero  back  scatter  is  obtained  can  be  found  using  equations  (81)  and  (82).  First,  equations  (81) 
and  (82)  are  multiplied  by  the  electric  and  magnetic  dipole  modal  functions  (see  equations  (91) 
and  (92)).  Then  we  equate  equations  (81)  and  (82)  and  solve  for  the  reflection  coefficient  y 
which  results  in  the  following  expression 


r.(K-iy  „  (94) 

'  (k  +  l)2 

From  this  expression  it  can  be  shown  that  the  magnitude  of  y  is  always  less  than  or  equal  to  one 
which  means  there  is  a  passive  reactive  load  or  dissipative  load  that  can  reduce  the  backscattered 
field  to  zero.  In  this  case,  the  backscattered  field  depends  upon  both  k  and  a. 


5.2.6  Realization  and  Implementation  of  a  Zero  Backscatter  Huygens  Source 

In  this  section  we  discuss  the  realization  and  implementation  of  a  zero  backscatter 
Huygens  source.  To  realize  a  Huygens  source  we  must  select  of  design  an  antenna  that  is  capable 
of  exciting  and  electric  and  magnetic  dipole  mode.  There  are  a  large  number  of  electric  and 
magnetic  dipole  antennas  that  can  be  found  in  the  literature.  However,  there  are  far  fewer 
combined  electric  and  magnetic  dipole  antennas,  particularly  for  scattering.  Over  the  course  of 
this  project  we  have  evaluated  a  large  number  of  these  for  use  in  a  low  backscatter  antenna.  The 
most  promising  candidate  is  the  folded  multiarm  spherical  helix  (FMSH)  [27]  which  we  will 
discuss  in  the  following  section.  An  example  of  the  FMSH  is  shown  in  Figure  59  The  Folded 
Multi-arm  Spherical  Helix  (FMSH) 


Figure  59  The  Folded  Multi-arm  Spherical  Helix  (FMSH) 


5.2.6.1  Folded  Multiarm  Spherical  Helix  (FMSH) 

The  folded  multiarm  spherical  helix  (FMSH)  can  support  both  the  electric  dipole  and 
magnetic  dipole  mode  [27,28]  which  makes  it  very  attractive  for  use  as  a  Huygens  source. 
Additionally,  it  has  been  shown  in  [27]  that  the  FMSH  can  approach  the  Chu-Wheeler-Thal 
radiation  Q  limit  for  an  electric  dipole. 
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The  dominant  mode  of  the  FMSH  can  be  either  an  electric  dipole  (TM)  or  magnetic 
dipole  (TE)  depending  upon  how  the  antenna  is  excited.  If  the  top  and  bottom  sections  of  an  arm 
are  driven  180deg  out  of  phase  the  electric  dipole  mode  is  dominant  because  the  far  field 
contribution  from  the  Iv  current  components  cancel.  This  excitation  is  illustrated  in  Figure  60.  If 
the  top  and  bottom  sections  are  driven  in  phase  (same  potential)  the  magnetic  dipole  mode  is 
dominant  and  the  far  field  contribution  from  the  Ig  current  components  cancel  as  illustrated  in 
Figure  61. 

Electric  Dipole  Excitation 


Figure  60  Electric  Dipole  Excitation  of  the  Folded  Multi-arm  Spherical  Helix  (FMSH) 


Figure  61  Magnetic  Dipole  Excitation  of  the  Folded  Multi-arm  Spherical  Helix  (FMSH) 


When  it  comes  to  realizing  a  Huygens  source  with  the  FMSH  there  are  essentially  two 
options.  The  first  option  is  to  use  the  same  FMSH  to  excite  both  electric  (TM)  and  magnetic  (TE) 
dipole  modes.  As  illustrated  in  Figure  62  this  can  be  accomplished  by  using  two  feeds  where 
both  feeds  are  essentially  dipoles.  However,  this  configuration  can’t  be  used  to  achieve  zero 
backscatter  because  the  electric  and  magnetic  dipole  moments  are  parallel.  In  addition,  this 
arrangement  results  in  greater  excitation  of  unwanted  spherical  modes  than  option  2 
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Figure  62  The  Same  FMSH  Used  to  Excite  Both  Electric  and  Magnetic  Dipole  Modes 

For  the  second  option,  the  electric  and  magnetic  dipole  modes  are  excited  by  separate 
FMSHs  with  one  being  placed  inside  the  other  or  nested.  This  leads  to  a  configuration  where  the 
magnetic  dipole  FMSH  (M-FMSH)  can  be  inside  electric  dipole  FMSH  (E-FMSH)  or  vice  versa 
(see  Figure  63).  It  is  known  that  a  spherical  electric  dipole  stores  less  energy  inside  (lower  Q) 
than  the  corresponding  spherical  magnetic  dipole.  Therefore,  one  might  assume  the  optimal 
arrangement  might  be  the  M-FMSH  inside  the  E-FMSH.  We  will  discuss  this  in  more  detail  in 
section  5. 2. 6. 3. 


M-FMSH  inside  E-FMSH 


E-FMSH  inside  M-FMSH 


Figure  63  Two  Nested  FMSH  Used  to  Eexcited  Orthogonal  Electric  and  Magnetic  Dipole 

Moments 


5.2.6.2  Implementation  Issues 

In  realizing  a  Huygens  source  there  are  a  number  of  issues  that  make  implementation 
challenging  regardless  of  the  antenna  type  used  to  excite  the  magnetic  and  electric  dipole  modes. 
First  of  all,  mode  balance  is  extremely  important  for  achieving  zero  backscatter.  To  achieve  zero 
backscatter  we  require  equal  excitation  of  the  electric  and  magnetic  dipole  modes  which  can  be 
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difficult  because  most  antennas  cannot  strongly  excite  a  magnetic  dipole  mode.  That  is,  if  you 
have  two  antennas  are  the  same  electrical  size  and  one  excites  an  electric  dipole  mode  and  the 
other  a  magnetic  dipole  mode,  the  electric  dipole  mode  tends  to  be  much  stronger  (i.e.  a  small 
wire  dipole  and  small  wire  loop  is  classic  example).  Fortunately,  the  FMSH  can  be  tuned  to 
excite  a  magnetic  dipole  mode  with  similar  strength  as  an  electric  dipole  mode. 

Another  key  issue  is  the  excitation  of  unwanted  higher  order  modes.  While  in  theory  it 
may  be  possible  to  achieve  zero  backscatter  when  higher  order  modes  are  present,  from  a 
practical  point  of  view,  it  becomes  significantly  more  difficult.  Therefore,  it  is  important  to 
suppress  the  excitation  of  higher  order  modes.  This  is  not  typically  a  problem  for  electric  dipole 
antennas  but  it  is  for  magnetic  dipole  antennas  because  we  must  use  electric  charge  to  create  a 
magnetic  dipole  moment.  This  would  not  be  a  problem  if  magnetic  charge  existed  but, 
unfortunately  it  doesn’t.  A  classic  example  of  this  is  a  small  loop  of  electric  current  which 
excites  a  magnetic  dipole  moment.  However,  the  excitation  is  not  pure  in  that  there  is  always  an 
electric  dipole  moment  whose  strength  depends  on  the  electrical  size  of  the  loop. 

Another  important  issue  is  coupling  between  the  dipoles.  Ideally,  there  is  no  coupling 
between  an  electric  and  magnetic  dipole  when  they  are  co-located.  However,  the  presence  of 
higher  order  modes  can  lead  to  coupling.  Also,  the  feeds  used  to  excite  the  dipoles  can  produce 
coupling  and  there  can  also  be  coupling  between  the  feed  networks  of  each  dipole.  All  of  these 
make  achieving  the  proper  mode  balance  more  difficult.  In  the  case  of  feed  network  coupling,  it 
can  also  lead  to  additional  terms  in  the  scattered  field  that  must  be  compensated  in  some  manner. 

The  last  main  issue  is  the  relative  phasing  between  electric  and  magnetic  dipole  modes. 
In  order  to  achieve  zero  backscatter  we  must  maintain  a  90  degree  phase  shift  between  the 
dipoles.  Any  deviation  from  this  will  increase  the  backscatter. 


5.2.6.3  Huygens  Source  Realization  using  FMSH 

In  this  section  we  will  show  how  the  implementation  issues  discusses  in  section  5.2. 6.2 
can  be  addressed  using  the  FMSH  antenna  to  realize  a  Huygens  source.  Firstly,  achieving  the 
desired  mode  balance  is  essentially  reliant  on  controlling  the  resonant  frequency  of  both  dipoles 
and  matching.  For  the  FMSH,  there  are  a  number  of  variables  for  controlling  the  resonant 
frequency  which  are  electrical  size  (with  respect  to  the  free  space  wavelength),  number  of  arms, 
number  of  turns,  termination/connection  of  arms,  feed  dipole  length  (M-FMSH  only).  Therefore, 
achieving  the  proper  mode  balance  is  achievable  (see  Figure  65  and  Figure  67).  To  demonstrate 
this,  let  us  first  consider  the  case  were  the  M-FMSH  is  inside  the  E-FMSH  as  shown  in  Figure 
64.  For  this  case,  both  dipole  modes  have  been  tuned  to  resonate  around  an  electrical  size  of  ka  ~ 
0.395.  Figure  65  shows  how  the  power  radiated  is  distributed  amongst  the  mode  modes  excited 
by  the  Huygens  source.  These  results  were  obtained  numerically  using  the  integral  equation 
based  software  FEKO  [29],  From  Figure  65  it  is  evident  that  most  of  the  power  is  in  the  electric 
and  magnetic  dipole  mode  which  means  that  higher  order  modes  are  negligible  over  the 
frequency  range  shown  in  the  figure.  Furthermore,  it  is  evident  that  there  is  equal  excitation 
(equal  power)  of  both  modes  near  the  resonant  frequency.  Outside  resonant  region,  there  is  a 
large  modal  power  imbalance  between  electric  and  magnetic  dipole  modes  and  the  antenna 
almost  behaves  as  a  single  mode  antenna.  This  behavior  is  not  desirable  because  we  need  equal 
excitation  of  the  modes  to  achieve  zero  backscatter.  Therefore,  the  Huygens  source  in  Figure  64 
is  fairly  narrowband. 
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Figure  64  Magnetic  Dipole  FMSH  Inside  an  Electric  Dipole  FMSH 


ka 


Figure  65  Modal  Power  for  Magnetic  Dipole  FMSH  Inside  an  Electric  Dipole  FMSH 

Let  us  next  consider  the  opposite  case  where  the  E-FMSH  is  inside  the  M-FMSH  as  shown  in 
Figure  66.  Again,  in  the  resonance  region  we  have  equal  excitation  of  the  electric  and  dipole 
modes  as  shown  in  Figure  67.  Outside  the  resonant  region  there  is  a  small  modal  power 
imbalance  between  electric  and  magnetic  dipole  modes.  Therefore,  this  implementation  of  the 
Huygens  source  is  potentially  more  broadband  than  the  former.  However,  for  this 
implementation,  higher  order  modes  outside  resonant  region  are  not  as  negligible  as  the  former 
antenna.  The  main  issue  is  the  dipole  feed  for  the  M-FMSH. 
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ka 

Figure  67  Modal  Power  for  Electric  Dipole  FMSH  Inside  a  Magnetic  Dipole  FMSH 


When  it  comes  to  mitigating  the  excitation  of  higher  order  modes  in  the  FMSH,  the  electric 
dipole  version  functions  extremely  well  as  shown  in  Figure  68.  The  results  shown  in  Figure  68 
are  based  on  a  single  excitation  of  a  single  pair  of  arms.  However,  each  pair  of  arms  can  be 
excited  to  further  suppress  higher  order  modes  as  the  electrical  size  increases.  The  magnetic 
dipole  version  of  the  FMSH  functions  well  as  a  magnetic  dipole  but  the  dipole  feed  excites 
additional  modes  outside  the  resonant  region  as  shown  in  Figure  69.  This  is  an  inherent  problem 
with  the  magnetic  dipole  version  of  the  FMSH  but,  to  some  degree,  this  can  be  addressed  by 
adjusting  length  of  feed  dipole  but  this  also  affects  the  resonant  frequency. 
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Figure  68  Modal  Power  for  FMSH  with  Electric  Dipole  Moment  Excited 
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Figure  69  Modal  Power  for  FMSH  with  Magnetic  Dipole  Moment  Excited 


We  know  demonstrate  that  it  is  possible  to  use  the  FMSH  to  realize  a  Huygens  source  with  zero 
backscatter.  In  Figure  70  the  total  scattered  field  in  the  yz  and  xy  planes  is  plotted  for  the 
Huygens  source  in  Figure  66  at  the  resonant  frequency.  As  you  can  see  from  Figure  70,  in  the 
backscatter  direction  (0=90  degrees)  the  scattered  field  has  been  greatly  reduced  but  it  is  not 
zero.  In  practice,  we  will  never  achieve  zero  backscatter  numerically  and  it  is  a  very  tedious 
process  that  involves  a  lot  of  fine  tuning  to  reduce  the  backscatter  to  a  very  low  level  but  it  is 
possible. 
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Figure  70  Scattering  Patterns  in  the  YZ-Plane  and  XZ-Plane  for  a  Huygens  Source 
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Acronym 

AFRL 

AMG 

CM 

CSWAP 

D 

dB 

EBG 

EM 

FEM 

FMSH 

FY 

GHz 

GRIN 

HFSS 

MHz 

NEMO 

Q 

PEC 

PMC 

PO 

Q 

STO 

TE 

TM 

TO 

UAS 

X 


List  of  Acronyms,  Abbreviations,  and  Symbols 

Description 

Air  Force  Research  Laboratory 

Artificial  Magnetic  Grounds 

Characteristic  Modes 

Cost,  Size,  Weight  and  Power 

Directivity 

Decibels 

Electronic  Bandgap 

Electromagnetic 

Finite  Element  Method 

Folded  Multi-arm  Spherical  Helix 

Fiscal  Year 

Gigahertz 

Graded  Index 

High  Frequency  Structural  Simulator 
Megahertz 

Numerical  Electromagnetic  Optimizer 
Ohms 

Perfect  Electric  Conductor 
Perfect  Magnetic  Conductor 
Physical  Optics 
Quality  Factor 

Source  Transformation  Optics 
Transverse  Electric 
Transverse  Magnetic 
Transformation  Optics 
Unmanned  Aircraft  System 
Wavelength 
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